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Abstract 

An operator theoretic approach to orthogonal rational functions 
on the unit circle with poles in its exterior is presented in this pa- 
per. This approach is based on the identification of a suitable ma- 
trix representation of the multiplication operator associated with the 
corresponding orthogonality measure. Two different alternatives are 
discussed, depending whether we use for the matrix representation 
the standard basis of orthogonal rational functions, or a new one with 
poles alternatively located in the exterior and the interior of the unit 
circle. The corresponding representations are linear fractional trans- 
formations with matrix coefficients acting respectively on Hessenberg 
and five-diagonal unitary matrices. 

In consequence, the orthogonality measure can be recovered from 
the spectral measure of an infinite unitary matrix depending uniquely 
on the poles and the parameters of the recurrence relation for the 
orthogonal rational functions. Besides, the zeros of the orthogonal and 
para-orthogonal rational functions are identified as the eigenvalues of 
matrix linear fractional transformations of finite Hessenberg and five- 
diagonal matrices. 
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As an application of this operator approach, we obtain new re- 
lations between the support of the orthogonality measure and the 
location of the poles and parameters of the recurrence relation, gener- 
alizing to the rational case known results for orthogonal polynomials 
on the unit circle. 

Finally, we extend these results to orthogonal polynomials on the 
real line with poles in the lower half plane. 

Keywords and phrases: orthogonal rational functions, unitary Hessenberg 
and band matrices, linear fractional transformations with operator coeffi- 
cients, pairs of operators. 
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1 Introduction 

The connection with Jacobi matrices has led to numerous applications of 
spectral techniques for self-adjoint operators in the theory of orthogonal poly- 
nomials on the real line. The direct extension of these ideas to the orthogonal 
polynomials on the unit circle yields a connection with unitary Hessenberg 
matrices (see [31 EEl EUJ [BO, [33]) which has provided some applications (see 
for instance [TH1 El UHl IH1 IBB] ) • Nevertheless, the authentic analogue of the 
Jacobi matrices for the unit circle is a class of unitary five-diagonal matrices 
which has been only recently discovered (see [361 H2])- This discovery has 
caused an explosion of applications of spectral methods for unitary operators 
in the theory of orthogonal polynomials on the unit circle, among which the 
numerous applications appearing in the monograph [30, 3TJ have been only 
the starting point. 

The orthogonal polynomials are a particular case of a more general kind 
of orthogonal functions with interest in many pure and applied sciences: 
the orthogonal rational functions with prescribed poles (see [10] and the 
references therein) . The natural generalization of the orthogonal polynomials 
on the real line and the unit circle requires the poles to be in the extended 
real line and in the exterior of the closed unit disk respectively. The first 
situation presents special complications, an indication of this being the fact 
that the poles can lie on the support of the orthogonality measure. Indeed, 
considered as orthogonal rational functions, the main difference between the 
orthogonal polynomials on the real line and the unit circle is not the location 
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of the support of the measure, but the relative location of the poles with 
respect to this support. Actually, the Cayley transform maps the orthogonal 
rational functions on the unit circle with poles in the exterior of the closed 
unit disk onto the orthogonal rational functions on the real line with poles 
in the lower half plane, so both of them can be thought as generalizations of 
the orthogonal polynomials on the unit circle. The purpose of the paper is 
to generalize to this kind of orthogonal rational functions the above referred 
spectral techniques for the orthogonal polynomials on the unit circle. 

An important ingredient in the theory of orthogonal rational functions 
are the linear fractional transformations z — > (a\z + CI2) (0.3,2 + 0,4) ~ 1 on the 
complex plane, where complex numbers. It is natural to expect the 

related spectral methods to have a close relationship with the operator ver- 
sion of such transformations, i.e., the maps T — > {A{T + A2)(A$T + A^) -1 
in the space of linear operators on a Hilbert space, where the coefficients Ai 
are now operators on the same Hilbert space. The theory of linear fractional 
transformations with operator coefficients goes back to the work [25] of M. 
G. Krein and Yu. L. Smuljan, motivated by the study of operators in spaces 
with an indefinite metric initiated by M. G. Krein in [23, 24J. As we will see, 
the matrices related to the rational analogue of the orthogonal polynomials 
on the unit circle are the result of applying a linear fractional transformation 
with matrix coefficients to the Hessenberg and five-diagonal unitary matrices 
associated with the polynomial case. 

This reason, and also a better understanding of the subsequent rational 
generalizations, motivates Section [2], which summarizes the basics on spec- 
tral methods for orthogonal polynomials on the unit circle and describes the 
main results needed about orthogonal rational functions on the unit circle 
with poles in the exterior of the closed unit disk. Section [3] introduces the 
operator linear fractional transformations of interest for such orthogonal ra- 
tional functions. The corresponding spectral theory is developed in Sections 
H] and [5j which are devoted to the approaches based on Hessenberg and five- 
diagonal matrices respectively. Section [6] presents some applications of the 
above spectral theory to the study of the relation between the support of 
the orthogonality measure and the poles and parameters of the recurrence 
relation for the orthogonal rational functions. Finally, the Appendix remarks 
the main analogies and differences with the spectral theory for orthogonal 
rational functions on the real line with poles lying on the lower half plane. 
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2 OP and ORF on the unit circle 



In what follows a measure on the unit circle T = {z G C : \z\ = 1} will be a 
probability Borel measure \i supported on a subset supp/x of T. Let /i be one 
of such measures and consider the Hilbert space of /z-square-integrable 
functions with inner product 

(f,g)n = J f(z)g(z)dfi(z), f,geLl. 

Unless we say the opposite we will suppose that supp/x is an infinite set. Then, 
{z n ) n >o is a linearly independent subset of L 2 ^ whose orthonormalization gives 
the orthogonal polynomials (OP) ((p n ) n >o with respect to /i. If we choose 
these polynomials with positive leading coefficient, they satisfy the recurrence 
relation 
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where (f^(z) = z n Tp n (l/z) and \a n \ < 1. This establishes a bijection between 
measures /ionT and sequences (a n ) n >i in the unit disk D = {z 6 C : \z\ < 1} . 

A central problem in the theory of OP on the unit circle is to find relations 
between the orthogonality measure /i and the sequence (a n ) n >i appearing 
in the recurrence relation for the OP. There are several approaches to this 
problem but these last years have seen a rapid and impressive development 
of new operator theory techniques (see [30j [311 [32] and references therein) 
based on the recent discovery of the analogue for the unit circle of the Jacobi 
matrix related to OP on the real line (see [T2l [36] ). 

The main tool for the operator theoretic approach to the OP on T is the 
unitary multiplication operator 

T»: Ll - Ll 

/(*) - */(*) 

It is known that the spectrum of coincides with supp/i and the eigenvalues 
of Tjj, which have geometric multiplicity 1, are the mass points of \l. The 
eigenvectors of a given eigenvalue A are spanned by the characteristic function 
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X{X} of the set {A}. Even more, if E is the spectral measure of T u then 
p(A) = (1, E(A)l) ^ for any Borel subset A of T. All these properties are 
true no matter whether suppp is finite or infinite. 

If {fn)n>o is a basis of L 2 ^, the matrix of T M with respect to (f n )n>o is the 
matrix M whose (i, j)-th element is My = (fi, T^fj)^. In other words, 



(zfo(z) zfx(z) 



fo(z) A(z) 



M. 



(2) 



Any matrix representation M of T M can be identified with the unitary oper- 
ator 

e^e 

x -+ Mx 

on the space t 2 of square-sumable complex sequences. This operator is uni- 
tarily equivalent to T M . Therefore, once we know the dependence of M on 
the parameters (a n ) n >i, this matrix permits us to recover the orthogonality 
measure /i starting from the recurrence relation of the OP. Regarding this 
problem, the utility of the matrix representation M depends on its simplicity 
as a function of the parameters (a n ) n >i. 

For instance, when the polynomials are dense in L^, the representation 
of Tjt with respect to the OP ((p n )n>o is the irreducible Hessenberg matrix 

(see ismsiEniEniES]) 
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Ji = (hi j) is called a Hessenberg matrix because hij = for i > j + 1, and 
the irreducibility means that ftj+ij 7^ for any j. Using the 2x2 symmetric 
unitary matrices 

On Pn 
Pn Of, 

the Hessenberg representation can be factorized as 



n > 1, 



(4) 



lim 

n 




e. 




(5) 
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where I and I n mean the identity matrix of order infinite and n respectively 
and the limit has to be understood in the strong sense. 

Apart from its complexity, the Hessenberg representation has the incon- 
venience of being valid only when the polynomials are dense in L 2 In the 
general case TC is a matrix representation of the restriction T M \ V: V — > V of 
T M to the T^-invariant subspace given by the closure V of the polynomials in 
L 2 . As a restriction of a unitary operator, T M \ V is isometric but not neces- 
sarily unitary. TL is a representation of T M iff any of the following equivalent 
conditions hold (see 115, 301): 



V = log y! <fil} m <* 



L n)n>l 



^ £ 2 <=>■ TC is unitary. 



We denote by m the Lebesgue measure on T. 

A way to avoid the problems of the Hessenberg representation is to use as 
a basis of L 2 U the Laurent OP (Xn)n>o that arise from the orthonormalization 
, . .), which are given by (see [121 EDI EH EE]) 



of (1, z, z -1 , z 2 , z 



X2n[Z) 



X2n+l{ z > = z V^n+AZ) 



n > 0. 



The corresponding representation of T M is the five-diagonal matrix (see 

W) 
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which, apart from being valid for any measure fi on T, is a band instead 
of a Hessenberg matrix. Also, it has a much simpler dependence on the 
parameters (a n ) n >i than in the Hessenberg case. Moreover, this five-diagonal 
representation has a much better factorization than the Hessenberg one since 
C = C a C e , where C Q and C e are the 2 x 2-block-diagonal symmetric unitary 
matrices 

e 3 _ e 2 



c n 



V 



e, 



V 



(7) 
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Alternatively, it is possible to orthonormalize (1 2 , ...)• This 

leads to the Laurent OP (Xn*)n>o where Xn*( z ) = X n (^/ Z )i i- e -; 

X2n*(z) = Z~ U ip 2n (z), X2n+l*(z) = ^Vl+li 2 ), U>0. 

The related representation of T M is simply the transposed matrix C T = C e C Q 
of C. 

The Hessenberg and five-diagonal matrices given in ([3]) and ([6]) represent 
a multiplication operator (on V or Lfy only when (a ra ) n >i lies on D. Never- 
theless, they are well defined matrices for any sequence (a n ) n >i in the closed 
unit disk D. Indeed, factorizations ([5]) and (171) show that, even in this case, 
the Hessenberg representation is isometric while the five-diagonal one is uni- 
tary. Furthermore, when some a n G T, these Hessenberg and five-diagonal 
matrices decompose as a direct sum of an n x n and an infinite matrix. 
This decomposition property is of interest when trying to make perturbative 
spectral analysis of such matrix representations. 

The Hessenberg and five-diagonal representations of T M also give a spec- 
tral interpretation for the zeros of the OP in terms of the parameters of the 
recurrence. This result comes from the relation between the OP and certain 
orthogonal truncations of the operator T M . The restriction T M \ V nt i of the 
multiplication operator to the subspace V n> i = spanjV, z l+l , . . . , 2;' +n_1 } 
has no sense since V n j, is not invariant under T M . To give sense to this kind 
of restriction we must multiply T M on the left by a projection on V n ,i- In 
particular, if P n y. — > L 2 ^ is the orthogonal projection on V nt i, the operator 
rpjJPnj) _ p n l T^ [- p n l j s ca lled the orthogonal truncation of T M on V n ,i- The 

key point is that, for any / G Z, the characteristic polynomial of T^ Pn ' l " > is, 
up to factors, the n-th OP <p n (see [30J). 

The first n OP (<^fc)fc=o are a Das i s °f 'Pn = V nt o and the related matrix of 
T^ n ' is the principal submatrix 7i n of H of order n. So, (p n is proportional to 
the characteristic polynomial of TC n , whose eigenvalues are therefore the zeros 
of <f n . Furthermore, for I = —[(n — l)/2], the first n Laurent OP (Xk)kZo 
are a basis of V n j and the corresponding matrix of ™ is the principal 
submatrix C n of C of order n. Hence, ip n is proportional to the characteristic 
polynomial of C n and, thus, the zeros of <p n are the eigenvalues of C n . 

Contrary to the full infinite matrix, Ti, n and C n are not unitary and depend 
only on the first n parameters oi, . . . , a n . However, factorizations (jSJ) and (J7j) 
show that if we change in these principal submatrices the last parameter 
a n G D by a complex number u G T, then we obtain a unitary matrix. The 
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corresponding characteristic polynomial is the result of performing n steps 
of recurrence (pQ), but substituting in the last one a n G © by u G T, i.e., it is 
a multiple of 

z<Pn-i(z) + Wn-lO)" 
Using (pQ), this polynomial can be alternatively written up to factors as 

<p n (z) +v<p* n (z), v= U _ , 

1 — a n w 

and the arbitrariness of u G T translates into a similar arbitrariness for 
v G T. These polynomials, called para-orthogonal polynomials (POP), were 
introduced for the first time in [21J. There it was proved that such POP have 
simple zeros lying on T, which play the role of nodes in the Szego quadrature 
formulas on T (the analogue of the Gaussian quadrature formulas on M), 
thus, providing finitely supported measures on T that *-weakly converge to 
the measure \i. Therefore, the nodes of the Szego quadrature formulas can 
be obtained as eigenvalues of Hessenberg or five-diagonal unitary matrices. 

Our aim is to generalize the above results to the orthogonal rational 
functions with poles outside of the support of the orthogonality measure. 
Two archetypical situations will be considered: measures on the unit circle T 
and measures on the extended real line R = 1U {00} ■ For convenience, the 
analysis will be done in a detailed way for measures on the unit circle, the 
discussion of the special features in the case of the real line being relegated to 
the Appendix. So, for the moment we will consider a measure /ionT and the 
corresponding orthogonal rational functions with poles arbitrary located in 
the exterior of the unit circle E = C \ D. We consider the extended complex 
plane C = CU{oo}to include for the poles the possibility of being located 
at 00. Indeed, the OP with respect to /i correspond to the special case of the 
orthogonal rational functions with all the poles at 00. 

An important transformation in C is z = l/~z, which leaves invariant any 
element of T and establishes a bijection between D and E. This transforma- 
tion induces the ^-involution f*(z) = f(z) in the set of complex functions, 
which defines an anti-unitary operator on L 2 ^ for any measure n on T. As a 
consequence, a sequence (/„) n >o of functions is a basis of iff (f n *) n >o is a 
basis too. Moreover, the *-involution on (T5]) gives 

(zf *(z) zf u (z) ...) = (/o*(z) fu(z) 



8 



which, taking into account that M is unitary, shows that the matrix of T M 
with respect to (f n *) n >o is the transposed M T of the matrix M associated 
with (f n )n>o- This relation holds when \x is finitely supported too, with the 
only difference that the basis of L 2 ^ are finite. 

Another essential ingredient in the theory of orthogonal rational functions 
on T are the Mobius transformations ( a defined for any a G © by 

r ( \ wl(z 




&J a (z) ' I vj*[z) = zw n Jz) = z — a. 



Up to factors in T, they are all the automorphisms of D. Indeed, ( a is a 
bijection of C onto C that leaves invariant T, D and E. The inverse trans- 
formation of ( a is ( a = C-a- It is a ls° remarkable that = l/( a . We 
distinguish the value a = that gives ( ao (z) = z. 

To get rational functions with fixed poles in E we introduce a sequence 
(«n)n>i in D. This sequence defines the finite Blaschke products (-B n )n>o 
given by 



Bq = 1, 



Bn Cai ' ' ' Ca n 5 Ti ^ 1 . 

Notice that B w = l/B n . The subspace 

C n = sp&n{B , B u ..., B n _i} = - 



(8) 



V 



consists of those rational functions whose poles, counted with multiplicity, lie 
on (ctk)k=i- We use the notation = span{£> n } n > = U n >iC n for the set 
of rational functions with poles lying on (a n ) n >i, counted with multiplicity, 
and C for the closure of in L 2 . 

If fi is a measure on T we can consider the rational functions ($„) ri >o that 
arise from the orthonormalization of (-B n ) n >o in L 2 ($„)„>o are called or- 
thogonal rational functions (ORF) with respect to fi associated with (a n ) n >i. 
When referring to ($ n )n>o we will call it in short a sequence of ORF on the 
unit circle. These functions satisfy a recurrence relation which, with an ap- 
propriate normalization of ($ n )n>o> has the form (see [TUl Theorem 4.1.3]) 

$o = l, 

<5>n\ _ (I b n \ f 2„Cn-l$n-l\ ^ ( 9 ) 

K J n ™n Un 1 K-l 1 ' - 
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where 



$*(«n-l)' ~" I 1 ifa n = 0, " \/ — „_j.(n„_| ) 1 |6„ 

and we use the notation 



Cn = Can, W„, = U7 an , W* n = CU* n , $* = Z X Z 2 • • • Z nJ B n $ n *- 

Notice that we do not follow the standard notation ( n = z n ( an and _B n = 
ziCai " " "^nCoin ( see f° r instance [TO])- In fact, concerning the matrix repre- 
sentations of the multiplication operator, it is more convenient to avoid the 
presence of the factors z n in recurrence Qj , something that we can get using 
the ORF (0„)„> o given by 

0o = 1, 

(j) n = z{z 2 -- -z n ^ n , n>l, 
and defining the superstar operation omitting the factors z n , that is, 



Then, flSJ) is equivalent to 

00 = 1, 



f I a n \ ( Cn-10n-l "\ ^ •, ( 10 ) 



with 



_ n (a n _i) / w n (a n ) 

a n — ~7 r — Zl%2 ' ' ' Z n n , 



0*(a n -i) " " " V w n-i(oi n -i) 1 - \a n \ 



2 ' 



In the polynomial case, corresponding to a n = for all n, ( ITUl) gives exactly 
(JTJ). As in the polynomial situation, the parameters (a n ) n >i of ( TTDT) lie on 
D. A Favard-type theorem also holds (see [10J Theorem 8.1.4]): given a 
sequence (a n ) n >i in D, the functions (0 n )n>o defined by recurrence (fTUjl are 
orthonormal with respect to some measure on T. This measure is unique 
when the infinite Blaschke product B(z) = Ti™=i Cn{z) diverges to zero for 
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z G D, i.e., when 52nLi(l — l tt n|) — oo. This condition means that the 
sequence (a„) n >i can not approach to T very quickly. 

Notice that, given a measure fi on T and a sequence (a n ) n >i in D, the 
parameters (a n ) n >i are uniquely defined. To see this, suppose that (0 n )n>o is 
another sequence of ORF satisfying a recurrence like (fTUj) . but with param- 
eters (a n ) n >i instead of (a n ) n >i. Then, with e n G T and e = 1. 
Hence, comparing the recurrences for (4> n )n>o and (0 n )n>o gives 

1 fe n 0\/l a„\ /e n _i \ = 1 / 1 a„\ 
^1 - |a„| 2 V ? n/ 1 / V e n-i) y/l - \a n \ 2 \&n I J' 

Taking determinants in both sides of the above equality we obtain | 

In consequence e n = e n _i for n > 1, which yields e n = e = 1. Therefore, 

a n = a n and n = <f> n . 

The above results show that any sequence a = (a n ) n >i in D defines a 
surjective application 

S a : — ► D°° 

fi — > a = (a n ) n >i 

between the set ^3 of probability measures on T and the set D°° of sequences 
in D. Furthermore, S a is a bijection when Y^n=i(^ ~~ \ a n\) — °o. The study of 
the application is one of the main interests to find a matrix representation 
of the multiplication operator T M with a simple dependence on the parameters 
a = (a n ) n >i and a = (a n ) n >i. Indeed, in the polynomial case, corresponding 
to et = 0, the five-diagonal representation C = C(a) of T u given in (jSJ) has 
revealed to be a powerful tool in the study of So. 

To find such a matrix representation, it is convenient to write recurrence 
(HOP in a different way. For any a G D we can define the positive number 

rj a = w a (a) l/2 = y/l - \a\ 2 . 

Denoting i] n = 7] an and introducing the parameters 

n i — m + Vn-l 

Pn=V 1_a n : Pn= Pn, P, 

Vn 

(HUD yields 

W n-\§n-\ = Pn^n4>n ~ On^n-l^-i, 
W n (j)* n = a n w n (p n + p~W n -x4>n-U 

This way of writing ([TO]) will be useful later. 



Vn-l 



Pn, (11) 



n > 1. (12) 
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3 Operator Mobius transformations 



As we will see, the operator version of the scalar Mobius transformations ( a 
appears in a natural way in the spectral theory of ORF on the unit circle. 
Analogously to the scalar case, such operator Mobius transformations are a 
particular case of the linear fractional transformations with operator coeffi- 
cients introduced by M. G. Krein in [231 121] for the study of spaces with an 
indefinite metric. A detailed study of these operator Mobius transformations 
in the general context of linear fractional transformations can be found, for 
instance, in the original paper of M. G. Krein and Yu L. Smuljan [25] or 
in the most recent survey of T. Ya Azizov and I. S. Iokhvidov [I] and the 
references therein. We will introduce the operator Mobius transformations 
summarizing the main properties of interest for us. 

Before doing this, we will fix some notations and conventions for linear 
operators. In what follows (H, (•,•)) means a separable Hilbert space. Given 
a linear operator T on H, denotes its adjoint, <j(T) its spectrum and cr p (T) 
its point spectrum. As usual, we omit the identity operator 1 on if so we 
use the same symbol z for the complex number z G C and for the operator 
zl, the meaning being clear from the context in any case. In general, we will 
deal with the Banach space (B#, || ■ ||) of everywhere defined bounded linear 
operators on H. 

In particular, M^n and can be identified with the sets of n x n com- 
plex matrices and infinite bounded complex matrices respectively. In this 
identification we associate any bounded square matrix M with the opera- 
tor x — > Mx, where x is a column vector of C n or I 2 . However, we could 
also consider the operator x — *> xM, where x is a row vector of C n or C 2 . 
Both operators have the same spectrum, although their eigenvalues can be 
different in the case of I 2 . Nevertheless, we will normally work with nor- 
mal or finite-dimensional matrices, for which the eigenvalues are the same 
in both situations. However, even in these cases, the eigenvectors are in 
general different. So, we will distinguish between right eigenvectors (or just 
eigenvectors) for x — > Mx and left eigenvectors for x — > xM. That is, right 
eigenvectors are the standard ones while left eigenvectors are the transposed 
of the eigenvectors of M T (in particular, when M is normal, right eigenvec- 
tors are the adjoints of left eigenvectors). In the subsequent discussions, this 
convention often permits us to avoid the T superindex, something convenient 
because many indices appear later. 

The operator Mobius transformations on H are linear fractional trans- 
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formations with operator coefficients that transform bijectively the unit ball 
Dff = {T G B>h : \\T\\ < 1} of B# onto itself. The role of the complex 
parameter a G D of ( a is played by an operator A G D#, so that 



Va 



V 7 ! - AAt 



w A (T) 


= 1 - TA\ 


w\(T) 


= T-A, 


w A (T) 


= 1 + A+T, 


zb* A {T) 


= T + A 



defines a positive operator with bounded inverse. Therefore, for any operator 
T in the closed unit ball B# = {T G B# : ||T|| < 1} we can define the 
operators ( A (T),( A (T) G M H by 

U{T) = i lA w A {T)- 1 w\{T)r A ], 

U(T) = r]- A l w* A {T) w A (T)- l VAU 

As in the scalar case, r\ A = w A (A) l l 2 . As we will see, the spectral theory of 
ORF is related to transformations ( A , ( A with A normal, so that r) A t = rj A in 
such a case. 

The transformations ( A and ( A are the operator analogs of the scalar 
Mobius transformations ( a and ( a respectively. The factors rj A , r) A x disappear 
in the scalar case due to the commutativity. Nevertheless, these factors are 
necessary for these operator transformations to keep similar properties to 
the scalar ones. Actually, ( A and ( A map B# on ID)#, as follows from the 
identities 

w A (T) 7^(1 - Ca(T) a(T)t) n1 x w A (T)t = 1 - TTt, 
w A (Ty V -l(l - UtYUT))vJ^a(T) = 1 - TtT. 13 

Besides, for any S,T G D#, a direct calculation shows that S = (a(T) iff 
T = Ca('S')) so ( a and Ca are mutually inverse transformations that map H) H 
onto itself. Furthermore, ([TBI also proves that ( A and Q A leave invariant B# 
and IV = {T G B# : ||T|| = 1}, mapping onto itself the set of isometries as 
well as the set of unitary operators on H . Indeed, as it was proven in [23] , up 
to unitary left and right factors, these operator Mobius transformations are 
the only linear fractional transformations with operator coefficients mapping 
bijectively B# onto itself. 

Using the relation rf A A = Arj At it is straightforward to verify the identities 

UtY = Gt(Tt), c ~ A(T) t = ^ t(T t )? (14) 
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which imply that (a(T) = U(T^ = £_a(T) as in the scalar case. Notice 
that the equalities (a = C-a and (a = (-a provide alternative expressions 
for ( A and ( A - 

Some formulas for the operator Mobius transformations will be of inter- 
est. From the relations (r] A ) n A = A{rj A] ) n for n — 0, 1, 2, . . ., and using the 
functional calculus for self-adjoint operators, we find that 

T] A A = Ar] A u 

Thus, if we define 

T A = VA lTr lAt 
for any linear operator T on H, then, for all T G B#, 

Ca(T a ) = w A {T)~ l w\(T), ( A (T) = w* A (T A ) w A (T A )- 1 . (15) 

This, together with the immediate identity 

w* A (T) - w A (T) S = Twa(S) - w* A (S), (16) 

yields 

w A {T) (C A (T A ) - S A ) = {T- US)) ^a(Sa) (17) 
for all T, S G B#. Substituting S by (a(S) in (fT7l) gives 

T — S = w A {T) r A \UT) ' US)) rij m- A (S), (18) 
where we have used that 

^a(USa)) = ^a(CUSa)) = W^ A {S)~ l r]%. 

If we take A = a and S — z with aGD and z£D, ffT8l) becomes 



z — T — (C«(*) - Ca(T)) «7«(T). (19) 

In particular, choosing T = A with AG©, 

C«(*) - Ca(A) = ^ ia) n , (* ~ A)- (20) 

Notice that (TT9"j) and ( 120|) actually hold for any z, A G C \ {«}. 
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4 ORF and Hessenberg matrices 



In this section we will prove that the orthogonality measure of a sequence of 
ORF, as well as the zeros of the ORF, have a spectral interpretation in terms 
of Hessenberg matrices. Our first aim is to find the matrix representation 
of a unitary multiplication operator with respect to a basis of ORF. Before 
stating the result, let us see which kind of matrix representation we can 
expect. Let v be a measure on T and (<p n ) n >o the corresponding OP with 
positive leading coefficient. Given a G D, the functions cj) n (z) = (p n (( a (z)) 
define a sequence (0 n ) n >o of ORF with fixed poles at a. The corresponding 
orthogonality measure is \i = u a , where u a (A) = u(( a (A)) for any Borel 
subset A of T. It is straightforward to see that recurrence (jTJ) for (y?„) n >o is 
rewritten in terms of (4> n )n>o as recurrence (fTUj) with the same parameters 
<*> = (fln)n>i, i-e., a = So(u) = S a ([i). The matrix representation of the 
isometric operator T v \ V with respect to (<^ n )n>o is a Hessenberg matrix 
Ti = 7~C{a) with the form ([3]). Therefore, the matrix of \ C with respect 
to the ORF (0 n ) n >o is 

(<<M*),^(*)U~=0 = «V<W,Ca(*)Vi(«)>v)S=0 = C«(w)- 
The following theorem is a natural generalization of this particular situation. 

Theorem 4.1. Let a = (a n ) n >i be compactly included in D ; /x a measure 
on T and a = (a n ) n >i = S a (n). Then, C is T^-invariant and the matrix of 
the isometric operator T M \ C with respect to the corresponding ORF (<ft n ) n >o 
is V = Ca(W), where H = 7i(a) is given in and 

\ 
J 

The isometric matrix V represents the full operator T M iff any of the following 
equivalent conditions is fulfilled: 

£ = <^> V = Ll <^> log \i ^ L x m ^ a £ f <^> V is unitary. 

Proof. \\A\\ < 1 because a is compactly included in D, thus £4 maps onto 
theirselves the sets of infinite isometric and unitary matrices. Therefore, 



A = A{ol) 



a 1 



Ci2 
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taking into account that 7i is isometric, CaC^O is a we U defined isometric 
matrix too. 

The starting point to prove the theorem is recurrence ([TO]) written as ()12p . 
The second relation in (fT2l) yields 



n-l 



tu n 0* = a n tu„0 n + 2J P n P„_i • ■ • P k+1 ak^k(fik, n > 1, 



(21) 



fc=0 



where we set a = 1. This identity, together with the first relation in ffT2"j) . 
gives 



h 



k,n 



k=0 

- a n+lPnPn-l * * * Pfe+l^fc if * < n, 

-a„ + ia n if = n, 

Pn+i if k = n + 1, 

if fc > n + 1. 



(22) 



w *a(. z ) ~ w a{ z ) Ti) = 0- 



If we define the matrix it = (hij), equality ( 1221) can be written as 

(M z ) M z ) ■■ 

Using (fTTj) we find that the Hessenberg matrix 

/-ai -pr a 2 -pip2 a 3 -pip 2 p3 a i -A 

— ai<2 2 —Hip^a^ —CI1P2P3CI4; 
TC — — 02^3 ~ci-2Pz a "i 

— a 3 a4 

v y 

can be related to the isometric Hessenberg matrix 7i given in ([3]) by 

H = Va 1 ^ Va = Ti-A, 



(23) 



(24) 



(25) 



where we have used that rj^ = tja since A is diagonal, so normal. From this 
relation, (fT51) and (IT6l we see that (1231) is equivalent to 



>o(2) M z ) 



z-(a(H)) =0. 



(26) 



This equality implies that C is invariant under T^, so the restriction T M f" jC is 
well defined and (aO~Q is its matrix representation with respect to (0„)„>o- 
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\ £ is an isometry because it is the restriction of a unitary operator, 
which agrees with the fact that CaC^O is isometric. Also, Ca(W) and W are 
unitary at the same time, that is, when a ^ £ 2 . Besides, T M \ £ is unitary 
iff T^£ = £. This implies that T™£ = £ for any n G Z, so {z n } n& i C £. 
Hence £ = L 2 ^ because span{z n } ng ^ is dense in L^. Conversely, if £ = L^, 
then T M f £ = is unitary. Therefore, Ca(W) is unitary iff the ORF (0 n ) n >o 
are a basis of L 2 i.e., iff (a(T~Q represents the full operator T M . Finally, it 
is known that the condition Y^Lifi ~ \ a n\) — 00 ? which is satisfied for a 
compactly included in D, ensures that £ = V (see [TUl Theorem 7.2.2]) and 
so it implies the equivalence between £ = Lj;, V = L 2 ^ and log// ^ L} m (see 
[H Corollary 7.2.4]). 

□ 

Given a measure /1 on T, the parameters a = So,(fi) corresponding to the 
ORF (<f) n )n>o associated with a are in general different from the parameters 
a,(°) = So(fi) related to the OP (^p n )n>o- For instance, if a n = a for all n, the 
comments at the beginning of this section show that S a (fi a ) = Sq(h). Taking 
into account that S a is a bijection for a constant sequence a, we conclude 
that S a (fi) 7^ <S (fi). Therefore, the equivalence V = L 2 ^ <^> a £ 2 is not 
trivial in the general case since the known result in the polynomial situation 
is V = L\ & o(°) it 2 . 

Contrary to the polynomial case, the unitary matrix V of the multipli- 
cation operator with respect the ORF basis is not a Hessenberg matrix in 
general, but Ca(V) is a Hessenberg matrix, where Ca is an operator Mobius 
transformation constructed using all the poles of the related ORF. In the 
polynomial case A = 0, thus Ca(V) = V an d ^ becomes a Hessenberg ma- 
trix. 

As a consequence of Theorem 14.11 and the spectral properties of the uni- 
tary multiplication operator, we have the following spectral interpretation of 
the support of the orthogonality measure for ORF. 

Theorem 4.2. Let a. be a sequence compactly included in Hi, fi a measure 
on T such that log/i' ^ V- m and (4> n ) n >o the corresponding ORF. If 

V = 0(W), A = A(a), H = H(a), a = S a (jj,), 

and £ is the spectral measure of V , then \i = £\ % \. Besides, supp/z = o~(V) 
and the mass points of fi are the eigenvalues of V, which have geometric 
multiplicity 1. X is a mass point iff (4> n {^))n>o £ Given a mass point 
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A, the corresponding eigenvectors ofV are spanned by (0o(A) 0i(A) •••) 
anrfM{A}) = (E^,l0»(A)IT 1 - 

Proof. Under the hypothesis of the theorem, V is the matrix representation of 
the full operator T M with respect to (0 n )n>o- Hence, if E is the spectral mea- 
sure of Tft, then //(•) = (4> , E(-)4> ) ^ = Also, supp^ = cr(T M ) = cr(V) 
and the mass points of \x are the eigenvalues of T M , that is, the eigenvalues 
of V, which have therefore geometric multiplicity 1. If A is a mass point, we 
know that the characteristic function of {A} is a related eigenvector of 
T M , so, {cf) n , n = /i({A}) 4> n (X) is the (n + l)-th component of a corre- 
sponding eigenvector of V. This implies that (0 n (A))n>o G ^ 2 - Conversely, if 
A is an arbitrary complex number such that (0 n (A)) n >o £ ^ 2 ; relation fl26l) 
shows that ( 0q(A) 0i(A) • ■ ■ ) is a left eigenvector of V with eigenvalue A. 
Due to the unitarity of V, A G T and the above statement is equivalent to 
saying that ( 0o(A) (f>i (A) • • • ) is a (right) eigenvector of V with eigenvalue 
A. Therefore, A is a mass point of /i. Also, the identity 

oo oo 

M{A}) = (* {A} , X {X} ), = Y,(X{X},<Pn)Mn,X{X})» = $>({A}) 2 |0n(A)| 2 

n=0 n=0 

proves that /i({A}) = (E^o I0»(A)| 2 ) -1 - 

□ 

The fact that the representation V is not a Hessenberg matrix, but a 
Mobius transformation of a Hessenberg matrix, makes the rational case more 
complicated than the polynomial one. However, the Hessenberg structure can 
be kept if we formulate the spectral results in terms of pairs of operators. 

Remember that, given a Hilbert space H and two operators T, S G M H , 
the spectrum and point spectrum of the pair (T, S) are respectively the sets 

<t(T, S) = {A G C : T — XS has no inverse in B#}, 
<Xp(T, S 1 ) = {A G C : T — XS is not injective}. 

In the finite-dimensional case both sets coincide. The elements of cr p (T, S) are 
called eigenvalues of the pair, and the eigenvectors of (T, S) corresponding 
to an eigenvalue A are the elements x G H\ {0} such that (T — XS)x = 0. In 
these definitions it is assumed that, if A = oo, T — XS must be substituted 
by S. 
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With the above terminology, the isometric matrix V and the Hessenberg 
pair {^* a {'H a ),'Sj a {'Ha)) have the same spectrum and eigenvalues because 
w a I(Ha) ±x G B £2 when ol is compactly included in D. So, Theorem 14.21 
can be obviously rewritten substituting V by the pair (w A (Tl A ), ^aO^a))- 
Notice that, given an eigenvalue A, ( <j>o(X) 0i(A) ■ ■ ■ ) is a left eigenvector 
of the pair, i.e., 

(0 O (A) 0!(A) ■~)(w* A (H A )-\& A {H A )) = 0. 

Moreover, H, A = r) A 1 'H'r) A with r/^ 1 G due to the restrictions on ex. 
Therefore, Theorem 14.21 also holds substituting V by the Hessenberg pair 
(zu A (TC), zb A (7i)), but the left eigenvectors with eigenvalue A are spanned by 
(0o(A) 0i(A) ■■■)r lA \ 

4.1 Zeros of ORF and Hessenberg matrices 

Let fj, be a measure on T and a = So(fi). As we pointed out in Section [21 
the characteristic polynomial of the orthogonal truncation of on V n = 
spanjl, z, . . . ,z n_1 } is a multiple of the n-th OP related to /i. From this 
result, the relation between the zeros of the n-th OP and the eigenvalues of 
the principal submatrices of TC(a) follows. 

To obtain a similar result for the ORF associated with a sequence at we 
have to consider the operator multiplication by ( n in L 2 , i.e., 

UT.y.Ll^Ll 

f—*Cnf 

and the orthogonal truncation of Cn(^)t) on C n . This orthogonal truncation 
is defined by Cn(7)J (£n) = L n ( n (T^) \ C n , where the operator L n : -> is 
the orthogonal projection on C n . The following theorem is the starting point 
to identify the zeros of the ORF as the eigenvalues of some finite matrices 
related to A(ct) and TC(a). 

We remind that the n-th ORF <p n has the form 

4>n = — , Vn G V n+ i \ V n , 7T n = W 1 ---W n , 

and the zeros of <f> n , which are the zeros of the polynomial p n , lie on D (see 
PH Corollary 3.2.2]). 
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Theorem 4.3. Let ct be an arbitrary sequence in J}, fi a measure on T and 
4>n — Pn/^n the related n-th ORF. Then: 

1. If Z n is the set of zeros of 4> n , ( n (Z n ) is the set of eigenvalues of 
Cn(^/i) an d these eigenvalues have geometric multiplicity 1. 

2. Ifp n (z) oc ni^i^ — ^k), the characteristic polynomial of Cn(^) 

n 

n(^-c(A fc ))- 

k=l 

Proof, f E C n \ {0} is an eigenvector of Cn(^)i)^ with eigenvalue w iff 
(L n ( n — w)f = 0, that is, L n (( n — w)f = 0. This is equivalent to state 
that (Cn — w)f G C^ Cn+1 = span{0 n }, or, in other words, / oc n (C« ~~ w)~ x . 
Writing w = ( n (X) and using (1201) we find that this condition can be expressed 
as f(z) oc p n (z)(z — A) -1 / V n _i(z) with A G Z n . This proves item 1. 

Item 2 is equivalent to assert that the algebraic multiplicity m w of any 
eigenvalue w = Cn(A) of Cn(X^ is equal to the multiplicity of A as a root 
of p n . Since the geometric multiplicity of w is 1, m w > k iff there exists 
/ G £ n such that (L„Cn — = and (X n Cn — w)^ 1 f ^ 0. Analogously to 
the previous discussion, we find that these two conditions are equivalent to 
/ G span{0 n (C n -w)~ J }J =1 \span{0 n (C„-w)" J }5 ; ~ 1 1 , i.e., to / = p/Vn-i with 
p(z) G span{wi- 1 {z)p n {z){z - A)- J '}J =1 \ span^" 1 ^)^)^ - X)' 3 }^, 
as can be seen using (1201) again. Hence, by induction on we find that 
m w > k implies that the multiplicity of A as a root of p n is not less than k. 
Conversely, if the multiplicity of A as a root of p n is greater than or equal 
to k, f(z) = (f) n (z)(( n (z) - w)- k oc w k ~ 1 (z)p n (z)(z - X)~ k /-K n ~i(z) G C n and 
the above results ensure that (L n ( n — w) k f = and (L n ( n — ^ 0, so 

m w > k. 

□ 

The next step is to obtain a matrix representation of the orthogonal 
truncation Cn(^u) > so that we can S^ ve a ma trix version of the above 
theorem. In the following results the subscript „ on a matrix means the 
principal submatrix of order n of such a matrix. This notation will be used 
throughout the rest of the paper. 
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Theorem 4.4. Let ct be an arbitrary sequence in J}, fi a measure on T and 
(0n)n>o the related ORF. If A = A(ct) and H = Ti.(a) with a = S a ([i), the 
matrix of C, n {T^ Cn ' > with respect to (4>k)kZo is Cn(V^); where 

Proof. \\A n \\ < 1, thus Ca* maps Tc^ onto itself. Also, from (jSJ) we obtain 
the factorization 




In-2 \ ( In-1 

©n-i / \ ~a n 



all the factors being unitary except the last one which has norm 1, thus 
\\H n \\ = 1. Hence, V (n) = CuO^n) is well defined and ||V (n) || = 1. A similar 
reason shows that Cn(V^"- ) ) is a well defined matrix with norm 1. 
To prove the theorem, let us write the first n equations of (1231) as 



(<p (z) • • ■ <f> n -i(z) ) (w^(z) - W An (z) Hnj = b n W n {z)(j) n {z) , 

b n = P +(o o ••• o i)ec. 



(27) 



Then, identities ( |T5|) . (|T6|) and the equality H n = Vjl^-nVAn = (7~Ln)A n 
obtained from (1251) transform (1271) into 



(Mz) ■■■ <t>n-i{z)){z-V^) = c n w n {z)<t) n {z), c n eC n . (28) 
Using (fT5|) we get 

{fate) 0n-l(z)){(n(z)-(n(V {n) ))=d n( P n (z), d n G C™. 

Hence, taking into account that 



0/t if < n, 
if Jfe > n, 



we finally obtain 

(L n C„0 O ••• ^nCn0n-l ) = ( 00 0n-l ) Cn ( V ( " } ) , 

which proves that Cn(V^) is the matrix of L n ( n (T^) \ C n with respect to 

□ 
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Theorems 14.31 and 14.41 provide a spectral interpretation of the zeros of 
ORF in terms of operator Mobius transformations of Hessenberg matrices. 

Theorem 4.5. Let a be an arbitrary sequence in 3, fi a measure on T, 
(4>n)n>o the corresponding ORF, A = A(ot) and H = H(a) with a = S a (/j,). 
IfVM = CaMu), then: 

1. The zeros of <p n are the eigenvalues ofV^ n \ which have geometric mul- 
tiplicity 1. If X is a zero of <p n , the related left eigenvectors ofV^ are 
spanned by ( O (A) • • • <j) n -i(X) ). 

2. 4> n = — with p n proportional to the characteristic polynomial of . 

Proof. From Theorems 14.31 and I4.4[ the eigenvalues of Cn(V^) have geomet- 
ric multiplicity 1 and cr((' ri (V^ n ' ) )) = ( n (Z n ), where Z n are the zeros of n . 
Also, the characteristic polynomial of CnO^ ) is ilfc=i( z — Cn(Xk)), where 
Pn{z) = Uk=x( z ~ o-(UV^)) = CXV (n) )), so, bearing in mind that 
Cn is bijective, <r(V (n * ) ) = Z n . Furthermore, given an eigenvalue A of V^ n \ the 
corresponding eigenvalue ( n (X) of Cn(V (n - ) ) has the same geometric and alge- 
braic multiplicity. Therefore, n^il- 2 — ^fe) * s the characteristic polynomial 
of V (n) . Finally, if A is a zero of <p n , (J2SD shows that (0 O (A) ■ • ■ n -i(A) ) 
is a left eigenvector of V^' with eigenvalue A. 

□ 

From ( fl5l) and ( fl~6l) we know that 

w *aS z ) -^A n {z)iin = zw An (H n ) -w* An (H n ) = (z- V {n) )w An (H n ). 
This gives other alternatives to express p n as a determinant, like 

p n (z) oc det (w* An {z) - w An {z)H n ^j = det (zw An (H n ) - ro^C^n)) • 

The interest of the above expressions is that they show that p n can be calcu- 
lated as a determinant of a Hessenberg matrix. Furthermore, the last expres- 
sion provides a new spectral interpretation of the zeros of the ORF, related 
to the concept of the spectrum of a pair of operators. It shows that the zeros 
of 4> n are the eigenvalues of the Hessenberg pair (zu An (H n ), w An (7i n )). Also, 
according to Theorem 14.51 the left eigenvectors of (uJ An (H n ), zu An (H n )) cor- 
responding to an eigenvalue A are spanned by (</>o(A) • • • n _i(A) ). Tak- 
ing into account that 7i n = nj~^li. n ri An , the zeros of 4> n can be also under- 
stood as the eigenvalues of the Hessenberg pair (ro^ (W n ), zb An (H n )), the left 
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eigenvectors with eigenvalue A being spanned by ( 0o(A) ■ ■ ■ n _i(A) ) rj A . 
Indeed, 

p n (z) oc det (to* An (z) - w An (z)n n ) = det (zw An (H n ) - w An (H n )) . 

Apart from the sequence (4> n )n>o of ORF, another remarkable rational 
functions arise in the theory of ORF. They are the so called para-orthogonal 
rational functions (PORF), given by 

Q V n = <Pn + V<j>l VeT. (29) 

The PORF are the generalization of the POP to the rational case. Anal- 
ogously to the POP, the interest of the PORF Q v n relies on the fact that, 
contrary to the ORF <p n , it has n different zeros lying on T which, thus, play 
an important role in quadrature formulas and rational moment problems 
(see PHI Chapters 5 and 10]). These quadrature formulas associate with each 
PORF Q" n a measure [i v n supported on its zeros with a mass (Y^k=o l^fc(^) H^ 1 
at each zero A. Such quadrature formulas are exact in £ n _i£ n _i*. 

Due to the exactness of the quadrature formulas, the first n OP ((fik^Zo 
related to fi are also an orthonormal basis of the n-dimensional Hilbert space 
L 2 ^. The convergence properties of the quadrature formulas imply that, 
for any sequence (f n )n>i in T, the sequence (/i^ n )n>i of measures *-weakly 
converges to the orthogonality measure /ionT whenever 5^^L 1 (1— |a n |) = oo. 

A spectral interpretation can be also obtained for the zeros of the PORF. 
To understand this, let us write the PORF in an equivalent way. Using 
recurrence fllUp we find that 

Q'n = (! + a n v)e n ^-^- (Cn-10n-l + «C-l)> U = Ca n ( V ) ■ ( 30 ) 

Notice that the parameter u goes through the full unit circle as the parameter 
v does so. fl30l) shows that, like <f) n , Q v n is obtained from n steps of recurrence 
ffTUj) . but changing in the n-th step a n G by u = ( a „( v ) e ^- The analogous 
substitution in H n gives 

which is obviously a unitary matrix. Therefore, CinC^O * s unitary too be- 
cause Q An preserves the unitarity. 
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The following result provides a spectral interpretation of the zeros of 
the PORF Q v n in terms of the unitary Hessenberg matrix TC™, as well as a 
connection of such a matrix with the unitary multiplication operator T^v . It 
can be understood as a limit case of Theorems 14.41 and 14.51 

Theorem 4.6. Let a be an arbitrary sequence in D, fi a measure on T, 
{<p n )n>o the corresponding ORF, A = A(ct) and H = TC(a) with a = S a (ii). 
If Qn = 4>n + v4>* n ^ s the n-th PORF related to v G T and \x v n is the associated 
measure, then: 

1. The matrix of T^ with respect to ((fik)kZo ^ s 

v (B!U) = U(K). « = C«». 

2. The zeros of Q v n are the eigenvalues ofV^ n ' u \ If X is a zero of Q v n , the 
related eigenvectors of V^' 1 ^ are spanned by (</>q(A) ••• </>„_i(A)) . 

q v 

3. Q v n = — with q" n proportional to the characteristic polynomial ofV^ n ' u \ 

7T„ 



Proof. Using (J2TJ) in ([30]) we find that 

' ' 1 + a n v 



k=0 



in _ J -Upn-lPn-2 ' ' ' Pfc+l«fc if A; < 71 — 1, r , x 

^ ~ \ -«a n _x if A: = n-l, n ' 

This relation together with the first n — 1 equations of (123!) lead to the matrix 
identity 

(0oW ••• cj )n _ 1 (z))(wX l (z)-^A n (z)ni)=b n w n (z)Ql(z), 

b n = Pn _ (0 • • • 1 ) G C n . 

1 + a n t> 

where ^ = VaM^ = (H u n ) An . So, (USD and (USD give 

(0o(*) ■•■ 0„_iOz))(*-V (n:u) ) = c n ^ n (z)Q^), C6C" (31) 

= in Zy,, thus (l3Tj) implies that y( n; ") is the matrix of T^v with respect 
to (<ftk)k=o- ^he res ^ °f the statements are a consequence of this one and the 
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properties of the multiplication operators, similarly to the proof of Theorem 
14.21 Alternatively, they can be obtained directly from relation (1311) . the 
unitarity of y( n: ") and the fact that Q v n has n different zeros. 

□ 

Analogously to the comments after Theorem 14.51 if u = Ca n (v), 
ql{z) cc det (w* An (z) ~ ™A n W H u n ) = det {zw An {H u n ) - zu* An (K)) , 

which gives q v n as a determinant of a Hessenberg matrix too. The zeros of Q v n 
are the eigenvalues of the Hessenberg pair {vtr\ C^C), ^.a„(W^)), whose left 
eigenvectors with eigenvalue A are spanned by ( 4>o(\) ■ ■ ■ </> n -i(A) ) rfj^ . 

5 ORF and five-diagonal matrices 

Apart from the presence of operator Mobius transformations, there are some 
drawbacks in the spectral theory of ORF previously developed: the appear- 
ance of a Hessenberg matrix 7i instead of a band one, the complicated de- 
pendence of Ti, = 7Y(a) on the parameters a, and the fact that it represents 
the full multiplication operator T M only for certain measures \i on T. We will 
not be able to avoid the operator Mobius transformations because they are 
linked to the ORF, but the other problems can be overcome by choosing a 
different basis of ORF in Lf,. 

The key idea is to use, instead of the ORF (0„) n >o with poles in E, other 
ones whose poles are alternatively in E and D. For this pourpose we define 
the finite odd and even Blaschke products 

B o = B o = h 

B%, — ClC3 •• -C2n-1, B e n = C2C4 • ■ ■ C2n, 71 > 1. 

Consider the rational functions (Xn)n>o given by 

X 2n = B e n J* 2n , X 2n + i = B e n J 2n+1 , n>0. (32) 

Since = l/( n , the subspaces Ai n = span{x , ■ ■ ■ , Xn-i} are 

M 2n = B^_ u C 2n = spanjS^, B°, B\„ . . . , B*_ u , B°}, 
M 2n +i = B^C 2n+ i = span{5^, B{, B{^ . . . , B°, B^}, 
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i.e., M.2n and A4 2 n+i are the sets of rational functions whose poles, counted 
with multiplicity, lie on (a 1} a 2 , • • • , a 2n -2, d 2n -i) and (di, a 2 , • • • , d 2n -i, a 2 n) 
respectively. We will use the notation .Moo = span{x n } n > = U n >i.M n and 
M. for the closure oi M. 00 m_ L 2 . 

The orthonormality conditions <ft n -L£ n and (0 n , n )/i — 1 can be rewritten 
using (p* n as 0*J_C n £ n and (0*, 0*)^ = 1. Hence, the orthonormality of ((p n ) n >o 
is equivalent to X2n-LBn*(2n£>2n = M 2n , X2n+i^B c nif C 2n+1 = M 2n +i and 
(Xn,Xn)n = 1, i-e., to the orthonormality of {xn)n>o- The sequence (x„)„>o 
is therefore the result of orthonormalizing (B^, B°, B{^ B^ -Bf*, . . .) in L^. 
Hence, relation (132]) establishes a bijection between ORF associated with the 
sequences (a„) n >i and (ai, d 2 , a 3 , d 4 , . . .). We can consider also the ORF 
associated with the sequence (di, a 2 , d 3 , a 4 , . . .), i.e., the ORF that arise 
from the orthonormalization of (B%, B° u , Bf, B^, B%, . . .) in L\. This ORF 
are (xn*)n>o, which are related to (0 n )n>o by 

X2n* = B°^(j)2n, X2n+1* = ^+l*0 2n +l) Tl > 0. (33) 

As a conclusion, the two possibilities to generate ORF which alternate 
poles in E and D are related between them and, also, to the ORF with poles 
in E. The last ORF have been extensively studied, thus every known result 
for them can be easily translated to the first ones. This is an interesting result 
because the ORF with poles arbitrarily located in C\T = EUD are not so well 
known than the ORF with poles in E. However, surprisingly, we will use the 
above connection to obtain new results for ORF with poles in E using ORF 
with alternating poles in E and D. The basic idea is that the ORF (Xn)n>o 
provide new matrix tools for the analysis of questions concerning the ORF 
(4> n )n>o- The reason for this is the different nature of the recurrence satisfied 
by (xn)n>o, which, as we will see, is a five-term linear recurrence relation. 
We could think that recurrence ffTUj) for (4> n )n>o should be better because 
it is two-term, but the presence of 0* causes a non-linearity which is the 
origin of the difficulties to connect with linear operator theory. Alternatively, 
expanding </>* in the basis (4> n ) n >o we find the linear relation (j22l . but this 
is a recurrence with a number of terms that increases with n, giving rise to 
a representation of T M related to a Hessenberg instead of a band matrix. On 
the contrary, the five-term linear recurrence for the ORF (xn)n>o provides a 
matrix representation of T M in terms of five-diagonal matrices, as the following 
theorem states. 
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Theorem 5.1. Let a be a sequence compactly included in D, p a measure on 
T and a = S a (p). Then, the ORF (Xn)n>o associated with (ai, a 2 , «3, «4, • • •) 
are a basis of L 2 ^ and the matrix ofT^ with respect to (xn)n>o isU — Ca{C), 
where A = A(at) and C = C(a) is given in (0|). 



Proof. Since a. is compactly included in D, ||^4|| < 1 and, hence, £4 maps 
unitarity matrices into unitary matrices. Thus Ca{C) is a weu defined unitary 
matrix because C is unitary. 

Using (JT2D and (H2) we find that, for n > 1, 



(34) 



ro 2n-lX2n-l = B^_ u (p 2n ZU 2n 4>2n ~ 1 2n^ ! 2n-l<i ) 2n-l) = 

— B^P2 n (ptn+l zu 2n+l4>2n+l ~ O2n+1^2n02n) — 

— B^_ lj fl2n\0'2n-\' UJ 2n-\^>2n-\ + P2n-l G7 2n-202n-2) 

— P2nP2n+l tI7 2n+lX2n+l — P2n a 2n+l ro 2nX2n ~ 

— &2n-\ a 2n w 2n-\X2n-\ — P2n-l a 2n^2n-2X2n-2, 

W 2nX2n = B^_ u (a 2n W 2 n4>2n + p2n W 2n-l<P*2n-l) = 

— B^a 2 n(P2n+l W 2n+l4 ) 2n+l ~ ^2n+l^ '2n<t>*2 n ) + 

+ B ( j l _ u p 2n (a2n-l'^2n-l4>2n-l + P2n-l a7 2n-202n-2) = 
= a 2nP2n+l W 2n+lX2n+l ~ a 2n a 2n+l' cu 2nX2n + 

+ 0'2n-lP2n ZU '2n-lX2n-l + P2n-lP2n ZU 2n-2X2n-2, 



while 

w oXo = pi&iXi ~ a i^oXo- 

This is the five-term linear recurrence for (Xn)n>o? which can be written in 
the form 

0, 



(Xo(z) Xi(z) ••■){w* A (z)-w A (z)C 
where C is the five-diagonal matrix 



(35) 
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Using (TTTj) we find that C can be related to the unitary five-diagonal matrix 
C given in ^ by 

C = rf^C r) A = C A . (37) 

Bearing in mind (|T5|) and (fl6l) . this relation implies that (|35|) is equivalent 
to 

(Xo(*) --O^-UC)) =0, (38) 

which shows that Ai is invariant under T M and Ca(£) is the matrix represen- 
tation of I" .M with respect to (Xn)n>o- 

Similar arguments to those given in the proof of Theorem 14.11 prove that 
I" Ai is unitary iff Ai = I? . However, T M \ Ai is unitary whenever 
a is compactly included in D because in this case Ca(C) is unitary for any 
sequence a in D. Therefore, Ai = L 2 , i.e., the ORF (xn)n>o are a basis of 
L 2 , which implies that Ca(C) is a matrix of the full operator T M . 

□ 

Remark 5.2. We know that (Xn)n>o an d (Xn*)n>o are basis of L 2 at the same 
time, and the corresponding matrices of T M are related by transposition. 
Therefore, the previous theorem can be equivalently formulated saying that 
(Xn*) n >o is a basis of L 2 whenever a. is compactly included in D and, in this 
case, the related matrix of T M is U T . Notice that the second equalitiy in ( TT4T) 
implies that U T = (a(C T ) because A is diagonal. 

Theorem 15 . 1 1 states that, contrary to the case of the ORF (0„) n >o, (Xn)n>o 
and (Xn*)n>o are basis of 1? for any measure /x on T if (a n ) n >i is compactly 
included in O. Indeed, the completeness of (x«)n>o an d (Xn*)n>o i n L 2 holds 
even under a more general condition for a, as the next theorem shows. De- 
noting T* = {/* : / G J 7 } for any set T of complex functions, the problem 
is to find sufficient conditions for the equality Ai = L 2 ^ or, equivalently, 
M m = L\. 

Proposition 5.3. Let a be a sequence inB>, ft a measure on T and (Xn)n>o 
the ORF associated with (a 1; a 2 , a 3 , a^, . . .). If 

oo oo 

J^ 1 - |«2fc-l|) = ~ l a2fc D = 00 ' 

fe=l fc=l 

£/ien (x„)„> and (Xn*)n>o are 6ot/i fraszs o/i^. 
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Proof. Given an arbitrary sequence (3 = (/? n ) n >i i n C\T, let us use the 
notation £oo(/3) for the set of rational functions with poles in (3 — ($ n ) n >i, 
counted with multiplicity, i.e., 

^oo(^) = |J ^- , (39) 

where vj^z) = z. Also, let C{(3) be the closure of £oo(/3) in L 2 . Notice that 
C((3)* = C(f3) and M = C(a>i, d 2 , a 3 , d 4 , . . .). 
We will show that 

(*) Y^T=i( l ~ l a 2fc-l|) = oo =>> {z J } jeN C £(«!, d 2 , «3, «4, ■■■), 

(") ZlfcLil 1 - l«2fc|) = oo {z~ 3 }jeN C £(«i,d 2 ,a 3 ,d 4 , • • •)• 

This demonstrates the proposition since span-jVjjgz is dense in L 2 ^. 

Indeed, we only must prove (i) since it implies (ii). To see this, as- 
sume that (i) holds for any sequence (cn n ) n >i in D. Then, applying (i) to 
the sequence (a 2 , a%, a^, a 3 , . . .) we find that ^^(1 — |ct 2 fc I ) = oo ensures 
•{V}j e N C C(a 2 , &i, a 4 , d 3 , . . .) = £(&!, a 2 , d 3 , a 4 , . . .), which, applying the 
*-involution, becomes {z^jjeN C C(aci, d 2 , a 3 , d 4 , . . .). 

The conditions YlkLii^ ~ l a 2Ai-i|) = oo and — \ a 2k\) = oo are 

equivalent respectively to the divergence (to zero) in D of the Blaschke prod- 
ucts B° = YlT=i C2fe-i an d B e = YlT=i C2fc- Thus, all what we must prove 
is that the divergence of B° implies that z^ e £(a 1; d 2 , a 3 , d 4 , . . .) for any 
JEN. 

According to ([39]), {z j } jeN C £oo(ai, d 2 , ot 3 , d 4 , . . .) = Aloo if a 2 fc-i = 
for infinitely many values fceN. Hence, we only need to study the opposite 
case that, without loss of generality, we can suppose is a± = a 3 = ■ ■ ■ = 
&2s-i — and a 2 fc_i ^ for k > s. Then, {z, . . . ,z s } C Aloo and, since 
(/, f)n < H/lloo for an y / e Aloo, it suffices to prove that 

inf 11^-/(^)1^^0, Vj>s. 

To measure the L°°-distance between a polynomial and a subspace like Aln 
we can use the following result (see [TJ p. 243] or the more recent reference 

[na p. i5o]): 



mm 



z N + q(z) 



W-l) ■ • ■ iz - W T 



a = a max{|w fc |,l} 
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Therefore, if B° diverges, taking n > s, 



inf 

/ e M2n 



I s+m , „s+m-l , 

Z -+- d m -\Z -+- ■ 



+ ai z s+1 -f{z) 



inf 



y 2n+m—l 



n 



n—l. 



«2fc-l, 



n 

k=s+l 



"2fc-l 



and 



inf 

/ e M2n + 1 

E C 



s+m 



I ^ ^s+rrt-1 | 



inf 



+ 0\Z 

z 2n+m + 



s+1 



/(*) 



n 



fc=iv 



«2fe)n 



fe=s+l 



2: 



«2fc-l 



n 

k=s+l 



"2fc-l 



for any m G N. This result implies by induction on m that z s+rn G M. for 
any m G N. 

□ 

In the polynomial case A = so U = C becomes a five-diagonal matrix. 
However, in the general case U is not a band matrix but its Mobius transform 
Ca(M) = C is five- diagonal. This fact makes the rational case more compli- 
cated than the polynomial one but, as we will see later, the matrix IA can be 
also used in the rational case to transcribe certain properties of the measure 
\x into properties of the corresponding sequence a = S a ([i). 

As in the Hessenberg case, the previous theorem provides a spectral in- 
terpretation of the support of the measure [i. The arguments are similar to 
those given in the proof of Theorem 14. 2\ but now the restriction log // ^ L x m 
is not necessary. Hence, we obtain the following result. 

Theorem 5.4. Let a be a sequence compactly included in Hi, fi a measure 
on T, (0 n )„>o the corresponding ORF and (Xn)n>o the ORF associated with 
(ai,a 2 ,a 3 ,a 4 ,...). If 

U = ( A (C), A=A(a), C = C{a), a = S a (ji), 

and £ is the spectral measure oflA, then n = E^\. Besides, supp/i = o~(U) 
and the mass points of fi are the eigenvalues of U, which have geometric 
multiplicity 1. X is a mass point iff (Xn(A)) n >o £ £ 2 ■ Given a mass point 
X, the corresponding eigenvectors oflA are spanned by (xoW XiW ' ") 



and M{A}) = (E~olXn(A)p 



(E~oK(A)h 
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Analogously to the Hessenberg case, we can formulate the above spec- 
tral results in terms of pairs of band operators. Theorem 14.21 implies that 
supp/i = a(zb A (C), wj\{C)) and the mass points of \i are the eigenvalues of 
the five-diagonal pair (zu A (C), w^{C)). Also, given a mass point A, the left 

~ ~ —1/2 

eigenvectors of (zu A (C), w A {C)) are spanned by (xo(A) XiW '"jVa • 
Furthermore, the factorization C = C Q C e makes possible to formulate the 
above results using the tridiagonal pair 

(w* A (C), w A (C)) C\ = (C Q + Ml C\ + AC a ) 

instead of the five-diagonal pair. 

5.1 Zeros of ORF and five-diagonal matrices 

The previous results suggest that it should be possible a spectral interpre- 
tation of the zeros of ORF and PORF in terms of five-diagonal matrices. 
Similarly to the Hessenberg case, an important ingredient for this is the or- 
thogonal truncation of Cn(^)t) on M. n . Taking into account that M. n = Bf^C n 
with I = [(n — 1)/2], the following generalization of Theorem 14.31 is of interest 
to relate this truncation to the zeros of the ORF <p n . This generalization 
deals with the orthogonal truncation of Cn(^jj) on hC n , h e L 2 which is 
given by UT^ hC ^ = L h n Cn{T^) \ hC n , where L h n . L 2 ^ -> L\ is the orthogonal 
projection on hC n . 

Theorem 5.5. Let ol be an arbitrary sequence in B ; /i a measure on T and 
4>n — Pn/^n the related n-th ORF. Then, for any Borel function h: T — > T: 

1. If Z n is the set of zeros of <fi n , ( n (Z n ) is the set of eigenvalues of 
Cn(^) ; and these eigenvalues have geometric multiplicity 1. 

2. If 'pn(z) oc Ylk=i( z ~^k), the characteristic polynomial of Cn(^At) ^ s 

n 

II(z-C»(A*)). 

k=l 

Proof. The operator multiplication by h in L 2 ^ 

h{T,):Ll^Ll 
f—*hf 
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is unitary because h maps T on itself. When restricted in the following way 

/ — * hf 

it yields a isometric isomorphism V between C n and hC n . The orthogonal 
projection on hC n is = h{T^)L n h{T^)\ where L n is the orthogonal pro- 
jection on C n . Thus, the orthogonal truncations of CnC^/J on hC n and £ n are 
related by Cn(T l iY h£n ' = V CnC?At) , so they are unitarily equivalent. In 

consequence, they have the same eigenvalues and with the same geometric 
and algebraic multiplicity. Hence, the result follows from Theorem 14.31 

□ 

Taking h = Bf^, I = [(n — l)/2], in the previous theorem we find that 
it holds for the orthogonal truncation of Cn(^)t) on A4 n , i.e., CnC^t) = 
M n ( n (T^) \ M. n , where M n : L 2 ^ — > L 2 ^ is the orthogonal projection on M. n . 
To give a matrix version of this result we simply need a matrix representation 

Of Cn(T^ M -\ 

Theorem 5.6. Let a be an arbitrary sequence in Hi, fi a measure on T 
and (xn)n>o the ORF associated with (a\, &2, CK3, cLt, . ..). If A = A(ot) and 
C = C(a) with a = S a (fi), the matrix of ( n (T^)( Mn ^ with respect to (xfc)fc=o 
is (niU^), where 

w {n) = a(c n ), iiw (n) n = i. 

Proof. From the factorization C = C Q C e we find that C n = C on C en . Only one 
among the factors C on and C en is unitary, but the norm of the remaining 
factor is 1, so ||C n || = 1. Since C^ n leaves T C n invariant, = (,A n (Cn) is a 
well defined matrix with ||W^|| = 1. The same holds for ^ n (L{^ n '). 

To prove that ( n (U^) is the matrix representation of Cn(^)t) with respect 
to (xk)k=07 us consider first an odd n. Then, the first n equations of (I35I) 
can be written as 

(Xa(z) ■ ■ ■ Xn-i(z) ) (vj* An { z ) ~ ™An{z)Cn} = b n w n (z)xn{z) , b n e C n . 

( |3"7j) gives C n = rf^C n r\A n = (C n ) An . This, together with identities (TT51) and 
dUD, yields 

{Xo{z) ■■■ Xn-l(z)) (Z-U {n) ) = C n W n (z) X n(z), C n E C™. 
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Using (TT§j) we get 



and, taking into account that 



M nXk 



Xk ifk< n, 
if k > n, 



we finally obtain 



( M n ( n xo 



M n ( nXn -i) = (Xo ■■■ Xn-i)UU {n) )- 



This equality proves that C n (U^ n ') is the matrix of M n £ n (T M ) \ M. n with 
respect to {Xk)lZl- 

On the other hand, if n is even, we consider the orthogonal truncation of 



C n (T M ) on M n *, i.e., Cn(T^ M ^ = M n *Cn(T M ) r M w , where M^.Ll -> 



is the orthogonal projection on .M n *. Taking into account that C is unitary, 
the identity obtained by applying the *-involution on (1351) reads 



A similar reasoning starting from the first n equations of this equality proves 
that the matrix of Cn(T fi Y Mn *^ with respect to (xk*)kZo * s Cn(W* ), where 
lli n) = ( An (Cn)- Notice that flU} implies that W* (n) = because A is 

diagonal. 

Given a measure /i on T, for any n e N, the subspace A4 n only depends on 
the parameters «i, . . . , ot n -i of the sequence ck, so the same holds for the or- 
thogonal truncations CniT^)^"^ and Cn{T tl )^ Mn * ' . Therefore, concerning the 
spectral properties of these truncations we can suppose without loss of gener- 
ality that a. is compactly supported on D. Then, the matrix representations 
of with respect to (xk)k>o and (xk*)k>o are U and U T respectively. Hence, 
the representations of the orthogonal truncations C,n{T^ MnS) and Cn(^) 
with respect to (xfc)fc=o anc ^ (Xfc*)fc=o are the principal submatrices (( n (U)) n 
and (( n (U T )) n respectively. The fact that {C <n {lA T )) n = (( n (U)) n implies 
that, when n is even, the matrix of Cn{T^Y Mn ' with respect to (Xfc)fc=o i s 



Remark 5.7. The proof of the previous theorem also shows that the matrix of 
the orthogonal truncation Cn(7)i) ( ' M ' l * ) with respect to {Xk*)'kZl is Cn(W (n) ) T . 



( Xo*(z) Xi*{z) ■■■) {^* A ( Z ) - ^a(z) (C t ) a ) = 0. 



(40) 



□ 
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As a consequence of Theorems 15.51 and 15.61 we have the following spectral 
interpretation of the zeros of ORF in terms of Mobius transformations of 
five-diagonal matrices. 

Theorem 5.8. Let a be an arbitrary sequence in 3, fx a measure on T, 
(</>n)n>o the corresponding ORF and (Xn)n>o the ORF associated with the 
sequence (ai,a 2 , as, 6:4, .. .). Let A — A(at) and C = C(a) with a = S a {fi). 
IfUW =Ca»(C*), then: 

1. The zeros of <fi n are the eigenvalues ofW n >, which have geometric mul- 
tiplicity 1. If X is a zero of 4> n , the related left and right eigenvectors 
ofU^ 1 ' are spanned by X n {X) and Y n (X) T respectively, where 

X n = ffp-ij (xo " " " Xn-l ) j Y n — B°n] ( xo* " " " Xn-l* ) • 

2. (p n = — with p n proportional to the characteristic polynomial ofl4^ n '. 

Proof. The vectors X n (z) and Y n (z) are rational functions with the poles 
lying on E, so they can be evaluated at any zero A of <p n since A G D. 
Besides, X n (X),Y n (X) ^ because B%x%h = 4>2ki B kX2k~u = ^fc-i and C 
has its zeros in E. 

The proof of the theorem is similar to the case of Theorem I4.5[ the only 
difference concerning the identification of the eigenvectors. To obtain the 
left eigenvectors of let us consider the first n equations of ( 1351 for an 
arbitrary n G N. These equations can be written as 

(Xo(*) ■•• Xn-l{z))(w* An {z)-W An {z)Cn) = 

= b n W n (z)Xn{z) + d n W n+ i(z)Xn+l{z), 

where b n ,d n G C n are 

Pn-i S n _i ) odd n, 
■ 1 ) even n, 
odd n, 
1 ) even n. 



p+(0 ■■■ 

-Pn a n+l ( ■ 



PnPn+l (0 
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Writing Xn and Xn* in terms of n and (f)* n with the aid of (j3"2|) and (13"3"j) . and 
using the first equation of (IT2|) in the case of even n, (j4"Tl) reads 



(Xo(*) ■•■ Xn-i(^)) - w A n (z)C n ^ = p+w n (z)Bl(z)^ n (z)v n , 

■■■ p+_ x a n _ x ) oddn, ( 42 ) 
• • • 1 ) even n. 



Also, remember that (ITBl . (TIE}) and C n = r]AC n r) An imply that 

(* - WW) w An (C n ) = w* An (z) - w An (z) C n . (43) 

Therefore, if A is a zero of <p n , (11211 and (1431) show that X n (A) is a left 
eigenvector of with eigenvalue A. 

Proceeding in a similar way with the first n equations of fj40l) we find that 
Y n (X) is a left eigenvector of U^ T with eigenvalue A for any zero A of 4> n . 
Therefore, Y n (X) T is a right eigenvector of I4^ n >. 

□ 

For a unitary matrix, like V in the case a ^ £ 2 , y(™ ; ") or U, left and right 
eigenvectors are related by the t-operation. However, this is not the case 
of the matrices or U^ 1 '. Theorem 14.51 only gives information about the 
left eigenvectors of V^ n \ while Theorem 15.81 provides both, the left and right 
eigenvectors of Li^ n \ Apart from the simplest form of U^ n \ this is another 
advantage of using this matrix instead of for the spectral representation 
of the zeros of ORF. 

Concerning the form of the eigenvectors of lA^ n \ notice that the factors 
B[( n -i)/2] an d B° n / 2 ] m X n and Y n are necessary to avoid any problem when 
evaluating them on a point of D. However, if a zero A of <fi n does not coincide 
with any for k — 1, . . . , n—1, then we can take as left and right eigenvectors 
(Xo(A), • • • ,Xn-i(A) ) and ( Xo*0)> • • • , Xn-i*W f respectively. 

As in the Hessenberg case, there are other alternatives to express p n as a 
determinant. Indeed, from (TT5l) . (ITS]) and the identity C n = T]^ C n T] An , 

p n (z) oc det (w* An (z) - w An (z)C n ) = det (zw An (C n ) - w* An (C n )) . 

So p n can be calculated as a determinant of a five-diagonal matrix. Further- 
more, the last expression shows that the zeros of n are the eigenvalues of 
the five-diagonal pair {w* An {C n ), w An {C n )). The associated left eigenvectors 

with eigenvalue A are spanned by X n (X) rj^J 2 . 
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Besides, the factorization C n = C on C en permits us to express p n as a 
determinant of a tridiagonal matrix. If n is odd, C en is unitary, thus 

Pn (z) (X det (z(C\ n + A\C on ) - (Con + AnCl)) 

and the zeros of <j) n are the eigenvalues of the tridiagonal pair 

(C n + -AnCln, C\ n + A n C on ), 

which has the same left eigenvectors as (vj a (C n ), (C n )). On the contrary, 
Con is unitary for an even n. In this situation we can use the fact that US n ^ 
and IA^ T = (^(C^) have the same characteristic polynomial, and the left 
eigenvectors of one of them are the transposed of the right eigenvectors of 
the other one. Hence, 

Vn(z) oc det (w* An ( z ) ~ ^A n (z)Cl) = det (zw An (C T n ) ~ ™*A n ( C n)) 
and, bearing in mind that C^ = C en C on , 

p n (z) cx det (z(C ] on + A ] n C en ) - (C en + A n Cl n )) . 
So, the zeros of n are the eigenvalues of the tridiagonal pair 

(C e n + -A-nC^n, Cq„ + A n C en ) 

— 1/2 

and the left eigenvectors with eigenvalue A are spanned by Y n (X) t]a . 

The zeros of the PORF Q n have a spectral interpretation in terms of 
band matrices too. Such an interpretation has to do with the matrix repre- 
sentation of T n v with respect (Xfc)fc=o) which is an orthonormal basis of L 2 ^ 
due to the exactness of the quadrature formulas associated with fi n . Similar 
arguments to those appearing before Theorem 14.61 show that the zeros of the 
PORF should be related to the unitary matrix obtained from C n when 
substituting the parameter a n £ D by u G T. More precisely, we have the 
following result, which can be understood as a limit case of Theorems 15.61 
and EH 

Theorem 5.9. Let a be an arbitrary sequence in B>, \i a measure on T, 
(<Pn)n>o the corresponding ORF, A = A(a.) and C = C(a) with a = S a (/i). 
If Q n = 0n + V( Pn ^ s ^ e n-th PORF related to v G T and \i v n is the associated 
measure, then: 
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1. The matrix of T^ with respect to (xk)'k=o ^ 

2. The zeros of Q v n are the eigenvalues ofU^ n]U ^. If X is a zero of Q v n , the 
related eigenvectors oflA 1 ^'^ are spanned by (xo(^) • • • Xn-i(A) ) . 

3. Q v n = — with q v n proportional to the characteristic polynomial ofU^ n ' u '. 

Proof. As in the case of Theorem 14.61 it suffices to prove item 1 . For an odd 
n = 21 + 1, using (I30p in a similar computation to that of (I34p gives 

<- 2 Xn-2 = Pn-lPi^nBtMU— ^ Pn-l u ^n-lXn-l ~ 

1 ' 1 1 + a„v 



n 1 

ffl n-2 fl n-iro n -2Xn-2 — P„-2 a n-l tZ7 n-3Xn-3? 
Qn 

1 + a n v 

+ a n -2P n -ltt n -2Xn-2 + Pn-2Pn-l W n-3Xn-3, 



W n-lXn-l — 0>n-\PZ,' aj nB\ n l2\* -, , Z - _ a «-l M ^n-lXn-l + 



where w = Can! 11 )- These relations can be combined with the first n — 2 
equations of (135]) in the matrix identity 

(Xo(*) ••• Xn-lW)(«7X.W-^(«)C«) = 

with CI = V A l C nVA n = (CW Thus, using (TJ5D and (USD we find that 

(Xo(z) ■■■ Xn-i{z)){z-U^) =c n w n (z)Bf,(z)Q v n (z), c n EC n . 

Therefore, U^ n ' u ' is the matrix of T^v with respect to (xk)kZo because Q v n = 
in Lf. v . 

rn 

On the other hand, if n = 21 is even, proceeding in a similar way with 
f l30|) and (135]) we arrive at 



Xo*(z) ■ ■ • Xn-u(z) ) (z - U {n ^ T ) = c n w n {z)B° u {z) Q v n {z), 



so W( n; ") T is the matrix of T^v with respect to (Xk*)k=o- Consequently, the 
matrix of T^v with respect to (xk)kZo ^ s M^ n '' u - 

□ 
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The zeros of a PORF can be also interpreted as eigenvalues of a pair of 
band matrices. If u = Ca n ( v ), 

q v n {z) oc det (w* An (z) - w An (z) Q) = det (zw An {C u n ) - ^K)) 

gives q" n as a determinant of a five-diagonal matrix. The zeros of Q v n are the 
eigenvalues of the five-diagonal pair (-aJ An (C^), a74 n (C^)) and, given an eigen- 

value A, (xo(A) • • • Xn-i(A) ) r) A J spans the corresponding left eigenvec- 
tors subspace. 

We have also a factorization = C'^ n C^ n , where C" n and C"„ are the result 
of substituting a n by u in C on and C en respectively (this substitution actually 
takes place only in C on or C en , depending whether n is odd or even). C% n and 
C" n are both unitary, so 

p n (z) oc det {z(C:l + AiC u on ) - (CI + A*CM 

and the zeros of n are the eigenvalues of the tridiagonal pair 

(C u + A f Mt r Mt 4- A >[ C U \ 

which has the same left eigenvectors as {w An {C^ l ),'UJ_A n {C'^)). 

6 Applications 

In this section we will present some applications of the spectral theory pre- 
viously developed for the ORF on the unit circle. We will use the results 
involving five-diagonal matrices due to their advantages. The corresponding 
spectral theory associates with each sequence of ORF a five-diagonal unitary 
matrix C(a) depending on the parameters a = (a n )„>i of the recurrence re- 
lation, and a diagonal matrix A(ct) depending on the sequence cu = (a„) n >i 
which defines the poles a n . These band matrices keep all the information 
about the ORF since they generate the full sequence of ORF through the 
associated recurrence. The importance of these matrices is that they play 
the role of a simple short cut that connects directly the parameters a, ct to 
the ORF and the related orthogonality measure. 

An essential difference with the polynomial case is that the matrix directly 
related to the ORF and the orthogonality measure is not the five-diagonal 
one, but an operator Mobius transform of it, namely, U(a, ct) = CA(a)(C( a ))- 
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This introduces important difficulties when trying to apply the spectral the- 
ory to the rational case. However, in spite of these difficulties, the matrix 
tool U(a, a) becomes powerful enough to deal with hard problems even in 
the rational case. To understand the scope of the rational spectral theory, 
we will use it to solve some non trivial problems about the relation between 
the behavior of the sequences a, a and the properties of the corresponding 
orthogonality measure /j,(a,ot). The answers to these problems are known 
for OP, but the generalizations to ORF are new. 

The strategy will be to apply standard results of perturbation theory to 
the unitary operator on £ 2 defined by the matrix U (a, a). We will apply such 
perturbation results to the comparison of U(a, a) with another normal ma- 
trix, eventually with the form U(b, (3). A useful remark for these comparisons 
is that, for (3 compactly supported in D, U(b,(3) defines a unitary operator 
for any sequence b in D since, then, C(b) is unitary. However, U(b,{3) only 
represents a multiplication operator on T when b lies on D. When b n E T 
for some n we know that C(b) decomposes as a direct sum of an n x n and 
an infinite matrix. Taking into account that A(f3) is diagonal, a similar 
decomposition holds for U(b,f3). 

The results of operator theory that we will apply state that two operators 
T, S on H have some common spectral property provided that the pertur- 
bation T — S belongs to certain class of operators. We will deal with two 
kinds of perturbations: compact and trace class operators. Both are subsets 
of B# that are closed under sum, left and right product by any element of 
Mh and also under the f-operation, that is, they are hermitian ideals of B#. 
This fact is the key that permit us to use techniques of band matrices in the 
spectral theory of ORF, according to the following result. 

Proposition 6.1. Let 3 be a hermitian ideal ofM H . If A, B e B> H are normal 
and AB = BA, the condition A — B e J implies the equivalences 

T-Se3 & ( A (T) - Cb(S) e 3 & U(T) - Cb{S) e3, WT,Se B H . 

Proof. It suffices to prove the first equivalence because (a = (-a- Let 3 be 
a hermitian ideal of B#. The identities 

TiT 2 - S 1 S 2 = (T\ - SO T 2 + S 1 (T 2 - S 2 ), T' 1 - S' 1 = -T~\T - S)S- 1 
prove that 

Tj, Si e Mh, T i -S i e3 T 1 S 1 -T 2 S 2 e3, 
T~ l ,S~ l eM H , T-Se3 =>■ T- 1 -S- l e3. 
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Suppose now A, B G D^- normal such that AB = BA and A — B G 3. Then 
Va~Vb = BB* — AA^ G 3. The functional calculus for normal operators 
shows that r] A r] B = r) B r]A, so rj A - rj B = {va + VbYK^a ~ Vb) e 3 since 
(?M + 1 G Bff because //a and are positive with bounded inverse. 
If, besides, T,S G B# are such that T — 5 G 3, then ^(T) — vj b (S) = 
SB^ — TA 1 " G 3 and ro^(T) - w* B {S) = T - S + B - A G 3. In consequence, 
T - S G 3 Ca(T) - ( B (S) G J. Substituting in this result A, B by -A, -B 
and T, S by Ca(^), (b(S) respectively, we also find the opposite inclusion. 

□ 

Taking into account that A(a) is diagonal, the above result has the fol- 
lowing immediate consequence. 

Corollary 6.2. If 3 is a hermitian ideal o/B| and ex, (3 are sequences com- 
pactly included in D, the condition A(ac) — A((3) G 3 implies the equivalence 

C(a) - C(b) G 3 U(a, a) - U(b, (3) G 3 

for any sequences a, b in D. 

Besides, from the factorization C(a) = C (a) C e (a) we find that, for any 
ideal 3 of , 

C (a) - C {b), C e {a) - C e {b) E3 C(a) - C(b) G 3. 

In fact, many perturbation results for the five diagonal matrix C(a) are known 
due to the extensive use of this matrix during the last years for the spectral 
analysis of OP on T. 

The perturbation results that we will use are the invariance of the es- 
sential spectrum for normal operators related by a compact perturbation 
(Weyl's theorem: see [35] and [HI ES]), and the invariance of the absolutely 
continuous spectrum for unitary operators related by a trace class pertur- 
bation (Kato-Birman theorem: see [22| [7] and [8]). Given an operator T, 
its essential and absolutely continuous spectrum will be denoted o~ e {T) and 
cr ac (T) respectively. In the case of a normal operator, cx e (T) is constituted 
by the limit points of er(T) and the eigenvalues with infinite geometric mul- 
tiplicity. In particular, for any measure \i on T, cr e (T M ) is the set {supp/i}' of 
limit points of supp/^ and cr ac (T M ) is the support of the absolutely continuous 
part fi ac of /i. 
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There are several ways to characterize the compactness of an operator 
but, in the case of an operator represented by a band matrix, a very practi- 
cal characterization is available: compactness is equivalent to stating that all 
the diagonals converge to zero. If the matrix is not banded the convergence of 
the diagonals to zero is only a necessary condition for the compactness. The 
compactness can be also used to characterize certain properties of the essen- 
tial spectrum. For instance, given a unitary operator T, cr e (T) c {A 1; . . . , A„} 
iff (Ai — T) ■ ■ ■ (A„ — T) is compact (Krein's theorem: see [2] and [IS]). 

The trace class operators, i.e., the operators T such that y/T^T has finite 
trace, are more difficult to characterize, even if they are represented by a band 
matrix. Nevertheless, any infinite matrix (kij) that satisfies the condition 
Si j < 00 represents a trace class operator on £ 2 . 

Concerning the compactness and trace class character of U(a, ct)—U(b, f3), 
Corollary 16.21 implies that it is a consequence of the same property for 
A(ct) — A(f3) and C(a) —C{b). The diagonal matrix A(et) — A(f3) represents 
a compact operator iff lim n (a n — (3 n ) = 0, and is trace class iff ct — (3 G £ 2 . 
As for the compactness and trace class arguments for C(a) — C(b) in the 
applications that we will discuss, they follow the same lines as in |30j . 

As a first group of applications in the study of the dependence fj,(a, a), we 
will analyze the extreme behaviors corresponding to a sequence a converging 
to zero or (subsequently) to the unit circle. In what follows Lim n x n means 
the set of limit points of a sequence (x„) in C. 

Theorem 6.3. If a is compactly included in D ; then: 

1. lima n = =3- supp/i(a, a) = T. 

n 

2. lim|a„| = 1 =3- {supp/i(a, a)}' = Lim C n (—a n a n+ i). 

n n 

3. limsup|a n | = 1 =>• fi(a,ct) singular. 

n 

Proof. It is straightforward to check that, for any sequence ct in D, the ORF 
((f) n )n>o corresponding to the Lebesgue measure 

dm(^) = f = -1A z = e « tf G [0,2vr), 
2ir 2m z 

are given by 

0o = 1, 

4>n = Vn B n _ u U>1, 
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and satisfy recurrence fllOp with parameters a n = 0. Therefore, when a is 
compactly supported in D, the unitary matrix U(0, a.) represents the multi- 
plication operator T m . So, a(U(0, a)) = suppm = T. 

Now, suppose an arbitrary sequence a in D such that lim n a ra = 0. Then 
C(a) — C(0) is compact, thus U(a, a.) — U(0, a.) is compact too. Hence, Weyl's 
theorem implies {supp/i(a, a)}' = {suppm}' = T, that is, supp/i(a, a.) = T. 



If lim r , 
matrix 



\a n \ = 1, then C(a) —V(a) is compact, where T>(a) is the diagonal 
( -ax \ 



Via) 



-aia 2 



-a 2 a 3 



\ 



J 



Therefore, Proposition 16.11 implies that U(a,a) — (^ a ^(V(a)) is compact 
too. Notice that 



/Co(-ai) 



V 



(i(-aia 2 ) 



C 2 (-a 2 a 3 ) 



is diagonal and bounded, so it is normal and Weyl's theorem states that 
{supp/x(a, a.)}' = Lim„C„(-a n a n+ i). 

Finally, assume that limsup n \a n \ = 1. This means that there is a subse- 
quence (a n )neZ) IcN, such that lim ne ja„ = a 6 T. Without loss of gener- 



ality we can suppose ^2 



1 1/2 



< OO, SO that J2n£l(\ a n ~ a \ + Pn) < °0 



because p n < y/2\a n — a\. Let b be the sequence defined by 



a if n G I, 
a n if n (fi X. 



The condition J2 n <=:i(\ a ~ a «l + Pn) < 00 ensures that C a (a) — C {b) and 
C e (a) — C e (b) are trace class, so the same holds for U(a, a) — U(b, a). The 
Birman-Krein theorem states that supp/x ac (a, a) = <j ac (U(b, a.)), but the fact 
that b n e T for infinitely many values of n implies that U(b, a.) decomposes 
as a direct sum of finite matrices, so it has a pure point spectrum and, hence, 
it has no absolutely continuous part. Therefore, fi ac (a, ct) = 0. 

□ 
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We can also obtain general conditions for the invariance of {supp/z(a, at.)}' 
and supp^ ac (a, at). 

Theorem 6.4. If at is compactly included in D, then: 

1. lim(a n - /3 n ) = lim(a n - b n ) = => {supp,u(a, at)}' = {supp,u(&, (3)}'. 

n n 

2. 22 _ Pn\ + I a n - K\) < oo =>> supp/i ac (a, a) = supp/i ac (b, /3). 

n 

3. Ifb n = X n a n with X n G C, then: 

lim|A n | = limA ra+ iA n = 1 =>• {supp/i(a, at)} 1 = {suppyu(6, a)}', 

n n 

^(||A n | 2 - 1| + |A n+ iA n - 1|) < oo supp/i ac (a, a) = supp/^ ac (6, a). 



^. /3 n — a n+ N, b n — a n+ N 



f {supp/i(a, a)}' = {supp//(6,/3)}', 
\ supp/z ac (a, a) = supp/z ac (6,/3). 



Proof. First, notice that any of the hypothesis of the theorem ensure that 
f3 is compactly included in D when at satisfies the same property. Thus, 
the spectral theory that we have developed works for both sequences, at and 
(3. Concerning the notation, in what follows we write p n = \fl 



usually, and a n = a/1 — \b n \ 2 . 

To prove the first item, notice that the inequality 

\pn - o- n \ 2 < \p 2 n - a 2 n \ = \ \a n \ 2 - \b n \ 2 \ < 2\a n - b n \ 

implies that the conditions lim n (o; n — (3 n ) = lim n (a„ — b n ) = ensure the 
compactness of A(at)—A((3), C (a)—C (b) and C e (a)— C e (b). Inconsequence, 
U(a,at) —U(b,/3) is compact too and Weyl's theorem implies the equality 
{supp/i(a, at)}' = {suppyu(b, (3)}'. 

Suppose now J2 n (\a n — [3 n \ + \a n — b n \) < oo. If limsup n \a n \ = 1, then 
limsup n \b n \ = 1, so we conclude from Theorem 16.31 3 that supp// ac (a, a) = 
supp/i ac (b, (3) = 0. If, on the contrary, limsup n \a n \ < 1, then \a n \, \b n \ < r 
for some r < 1. Taking into account the inequality 

i | . IKI 2 - \b n \ 2 \ r 

\Pn - &n\ < ; < i = Qn ~ n , 

Pn + cr n VI - r 2 
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J2 n (\ a n- Pn\ + \a n — &n|) < °o implies that A(ex) — A((3), C (a)—C (b) and 
C e (o) — C e (b) are trace class. Thus, U(a, ex) — U(b, (3) is trace class too and, 
from the Birman-Krein theorem, suppp ac (a, ex) = supp/i ac (b, /3). 

1. We can write 
so that lim n \9 n+1 — 9 n \ = 0. 



Consider b n = X n a n with lim n |A 
A n = |A n |e i£ln with 6 n e [9 n -\ - ir,9. 
Define 

/ Ui 

u 2 

(' u 3 



lim n A„ + iA 

n 



, l/ n _i + 7T 
\ 



111 



V 



ii, 



••■/ 



The identity 



e„(a) 



■it. 



- e„(6) = 

a„u n (|A n | 2 

PnUn+lUn ' 



"I) 

On 



a n u n {u nJr iU n — |A n | 



together with lim n u n+ \u ri 



1 and \p n - a n \ 2 < \p\ 



°n I — 1 1 — I An | 1 1 On | , 

shows that UC (a) U—C {b) and Z7'C e (a) U'—C e (b) are compact. This implies 
the compactness of UC(a) W - C{b) = UC Q {a)U U ] C e {a)U ] - C {b)C e {b), 
which, bearing in mind Proposition 16.11 is equivalent to the compactness of 
UU(a, a)W - U(b, ex) = ( A(a) (UC(a)W) - ( A(a) (C(b)). Therefore, Weyl's 
Theorem implies that {supp/i(a, ex)}' = {supp/i(b, a)}'. 

When ^ n (||A n | 2 — 1| + |A„ + iA„ — 1|) < oo we have to consider again 
two possibilities. If limsup n |a n | = 1, necessarily limsup n |6 n | = 1 and 
supp/i ac (a, ex) = suppyti ac (b, a) = from Theorem 16.31 3. If limsup n |a„| < 1, 
then \a n \, \b n \ < r for some r < 1, so the relations 



\pn ~ er n | < 
I- 2 — 2 — 1 1 

\ U n+l U n ~ l \ 
\ u n+lU n — 1| 



- 16 l 2 l 

"n 



Pn + O'n 

' |A n +iA n 

\ U n+l U n 



< 



2V1 



I A, 



|A n +iA„|| < 2|A n+ iA n 

H 1, o _ 



\U n+ iU n + 1 



III 



n+l u n 



1| ^ |A ra +iA n — 1|, 



ensure that UC (a)U — C Q {b) and WC e (a)W — C e {b) are trace class. The 
Birman-Krein theorem then proves that supp/i ac (a, a) = supp/i ac (b, ex.) sim- 
ilarly to the previous case. 
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Finally, let b n = a n+ jy and (3 n = a n+ N for some iV G N. Consider the 
sequences a and a. given by 

1 if n < N, _ f if n < N, 

Ct — \ 

a n if n > AT, n \ a n if n > N. 

A(ct) — A(ct), C Q (a) — C Q (a) and C e (a) — C e (a) are finite rank, therefore 
U(a, a) — U(a, a) is compact and trace class. Besides, we have the decom- 
position U(a, a) = —I N @U(b, f3), so U(a, ct) and U(b, (3) have the same es- 
sential and absolutely continuous spectrum. As a consequence of these facts, 
the Weyl and Birman-Krein theorems give {suppyu(a, a)}' = {suppyLt(&, (3)}' 
and supp/i ac (a, a) = supp/x ac (6, (3). 

□ 

Combining the different results of the previous theorem we can obtain a 
more general one. 

Theorem 6.5. For any sequence ct compactly included in D; 

1. If \im(a n+N - (3 n ) = \im(\ n a n+N - b n ) = 0, lim |A n | = lim A„ + iA„ = I, 

n n n n 

then {supp/i(a, a)}' = {suppyu(6, /3)}'. 

2. If^2(\a n+N - (3 n \ + \\ n a n+N - b n \ + \ \\ n \ 2 - 1| + |A„+iA n - 1|) < oo ; 

n 

then supp/i ac (a, a) = supp/i ac (6, (3). 
A particular case of this theorem is worthwhile to be emphasized. 
Corollary 6.6. Let a e D, a e [0, 1], A e T and 

I\ a = {\e ie : \9\ < 2arcsina}. 

1. Iflima n = a, lim|a„| = a and lim n+1 = A, then 

n n n a n 

{supp//(a, a)}' = T \ C Q (r A ,„). 

2. If (\a n — ot\ + ||a n | — a| + " +1 — A j < oo ; t/ien 

supp// ac (a, a) = T \ Ca^a)- 
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Proof. Let us write a n = \a n \v n with v n G T. Notice that ex is compactly 
included in D because it is convergent in D. Therefore, we can apply Theorem 
16.51 to /i(a,a) and /x(6, (3) with (3 n = a, b n = X n a and A n = X n v n . Taking 
into account the relation 



I A n +lA n 



11 = 


A Vn+l 


< 


A an+l 


+ 


a n+l v n+l 















^ a n+l 


+ 







- 1 



< 2 



*n+l 



we find that {supp/i(a, a)}' = {suppyu(6, (3)}' under the assumptions of item 
1, and supp// ac (a, ex) = supp/z ac (6, /3) under the hypothesis of item 2. On 
the other hand, from the comments at the beginning of Section HJ we know 
that fi(b, (3) = u a , where u = fj,(b, 0) is the measure on T whose OP have 
parameters X n a and v a is defined by z/ Q (A) = u(( a (A)) for any Borel subset 
A of T. Therefore, {suppz/ a }' = (^({suppz/}'), supp(z/ Q ) ac = C a (suppz/ ac ) 
and the corollary follows from the well known result {suppz/}' = suppz/ ac = 

T\r\ a . 

□ 



If a = 0, Corollary 16.61 1 is a direct consequence of Theorem 16.31 1. while 
Corollary 16.61 2 can be derived from Szego's Theorem for OP on T: The- 
orem EH 2 implies that supp/i ac (a, a) = supp/i ac (a, (3) for f3 n = a when- 
ever ^2 n \ct n — a| < oo. n(a,f3) = u a , where now v = fi(a, 0), and the 
condition ^ n \a n \ < oo gives suppz/ ac = T because v is in the Szego class 
a = So(u) G £ 2 . Hence, supp/i ac (a, ex) = supp(z/ Q ) ac = T. In fact, this 
reasoning proves that the equality supp/i ac (a, ex) = T holds under the more 
general condition ^2 n (\a n — a| + |a n | 2 )<oo. 

Corollary 16. 6l of Theorem 16.51 can be understood also as an example of the 
following general result. It says that, when ex is convergent in D, the analysis 
of {suppyu(a, ex.)}' and supp/i ac (a, ex.) can be related to the much more known 
polynomial case, corresponding to ex = 0. 

Theorem 6.7. Let a G D. 

1. lima n = a =>- {supp/i(a, ex)}' = £ a ({supp//(a, 0)}'). 

n 

2. |a n - of] < oo =$> supp/i ac (a, ex) = C Q (suppyU ac (a, 0)). 



46 



Proof. Again ex is compactly included in © because it is convergent in D. So, 
if j3 n = a, Theorem 16.41 implies that {supp/i(a, ex)}' = {supp//(a, (3)}' when 
lim„a n = a, and suppp ac (a, ex) = suppp ac (a, /3) when \a n — a\ < oo. On 
the other hand, /x(a, /3) = u a with i/ = //(a, 0). As in the proof of Corollary 
16.61 the result follows from the relation between v and u a . 



The importance of the above theorem is due to the numerous known 
results for the relation between \x and a in the case of OP on T. Theorem 
16.71 permits us to translate some of these results to those ORF on T whose 
poles converge in E. For instance, Corollary 16.61 1 can be understood as the 
translation to this kind of ORF of a result for OP on T due to Barrios-Lopez 
(see [5]). This result was generalized later on in [26j as an improvement of 
a partial extension appearing in [T4]. The corresponding translation of this 
generalization to ORF states that Corollary 16.61 1 holds even if we substitute 
the condition lim n \a n \ = a by the more general one liminf n | 

All the above results provide only sufficient conditions on the sequences 
ex and a to ensure a certain property for the measure fi(a,ex). On the 
contrary, Krein's theorem permits us to characterize exactly those measures 
/x(a, a) with a fixed finite set {supp/i(a, ex.)}'. The characterization is in 
terms of the compactness of a matrix depending on a and ex. The fact that, 
contrary to the polynomial case, this matrix is not banded makes difficult to 
translate its compactness into equivalent conditions for the sequences a and 
ex. Nevertheless, in the case of {snppfi(a, ex)}' with at most two points we 
can find explicitly such equivalent conditions. 

Theorem 6.8. If ex is compactly included in D and A, Aj, A2 G T ; then: 
1. {supp /i(a, ex)}' = {A} iff 



□ 



limCn( 



«n«n+l) — A. 



n 



2. {supp/i(a,a)} / C {A 1; A 2 } iff 
lim n p„p„, + i = 0, 



( w n (\x) ^»-iQ 2 ) , , A 
hmp n — rfc„(A 2 ) -, r«„-i(Ai) = 

n \W n {a n ) W n -\\OL n -\) J 



lim (fc„(Ai)fc n (A 2 ) + (p-) 2 ^(Ai) U7;_ 1 (A 2 ) + 

n \ 
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where k n (z) = a n w* n (z) + a n+l w n (z). 

Proof. We are dealing only with measures \x with an infinite support on T, 
thus, supp/i has at least one limit point in T. Hence, from Krein's theorem, 
{supp/j}' = {A} iff A — T M is compact, i.e., iff A — U is compact. (fTTj) yields 

A - U = A - UC) = i£Va(A) (UA) - C) w x (C)- V (44) 

Bearing in mind that 77.4., 07.4(A) and 074(C) are bounded with bounded in- 
verse, the above expression shows that the compactness of A— U is equivalent 
to the compactness of Ca(A) —C. On the other hand, Ca(A) — C is compact iff 
\im n p n = and lim n (£ n (A) + a n a n+ i) = 0. However, the first of these condi- 
tions is a consequence of the second one because |Cn(A) + a n a n+ \\ > 1 — \a n \ 
since A G T. Also, taking into account (ITTj) . 

07 n (A) (Cn(A) + a n a n+ i) = (A - C„(-a n a n+ i)) ro n (-a n a n+ i). 

Therefore, lim n ((' n (A) + a„a„ +1 ) = iff lim„(A — Cn( — a n a n+ i)) = because 
2 > |ca n (A)|, |w n (— a n a n+ i) I > 1 — |ct n | and a is compactly supported in D. 

As for the case of two limit points, Krein's theorem implies that the 
inclusion {supp/i}' C {Ai, A 2 } is equivalent to the compactness of the matrix 
(Ai — W)(A 2 —IX). To express this condition as the compactness of a band 
matrix we use (jlij) for the factor A 2 — U, but for Ai — U we use the equality 

X — U — X- C- A (C) = VA^-A(Cr\CAW - C) w_ A (X) Va \ (45) 

obtained from (1T71) and the identity £4 = (-a- Then, similarly to the case of 
one limit point, we find that (Ai — W)(A 2 —IX) is compact iff the 9-diagonal 
matrix (O(Ai) — C) vjj^Xi) w^(^4) _1 w^(A 2 )(C_4(A 2 ) — C) is compact. This 
compactness condition can be equivalently formulated using a simpler band 
matrix obtained multiplying the above one on the left and the right by the 
unitary matrices C\ and C\ respectively. Taking into account the identity 
w* A {z) = zvoa{z)\ z G T, we find in this way that {supp/u}' C {Ai,A 2 } iff 
the five-diagonal matrix K{Xi)^wa{-A) 1 K{X2) is compact, where K(z) = 
vo* A {z) C\ — zua(z) Co- Now, it is just a matter of calculating the diagonals of 
K {Xi)^ 'w a(A)~ 1 K '(A2) to obtain the conditions given in the theorem. 

□ 

The implication lirn n £ n (— a n a n+ i) = A G T =>> {supp/j}' = {A} was in 
fact a consequence of Theorem 16.31 2. Krein's theorem adds the opposite 
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implication. Concerning the case of two limit points notice that, although 
the third condition is symmetric under the exchange of Ai and A2, the second 
one does not show explicitly such a symmetry. However, a detailed analysis 
of the second condition reveals that it is symmetric too. 

It seems that there is no simple way to generalize the arguments given 
in the proof of Theorem 16.81 to the case of more than two limit points. The 
reason is that, for n > 3, identities (jlij) and are not enough to reduce 
the compactness of (Ai —IX) • ■ ■ (A n —IX) to the compactness of a band matrix. 
So, contrary to the polynomial situation (see [16J and [131 [26]), the practical 
application of Krein's theorem to characterize in terms of the sequences a 
and ex those measures on T whose support has a finite set of more than two 
limit points remains as an open problem in the rational case. 



7 Appendix: ORF on the real line 

In what follows, a measure on the real line will be probability Borel measure 
/i supported on an infinite subset supp/i of R. When 00 is not a mass point 
of n we will refer to fi as a measure on R. Notice that we are considering all 
these measures as measures on R, no matter whether they have a mass point 
at 00 or not. This means that 00 G supp/i when 00 is a mass point of \x or 
when \i is a measure on R with unbounded standard support, so that supp/i 
is always closed in R. 

Analogously to the case of the unit circle, for any measure [i on the 
real line it is possible to consider ORF in L 2 ^ with poles in the lower half 
plane L = {z G C : Im(,z) < 0}. For this purpose we introduce for any 
a 6 D = 6 C : Im(z) > 0} the linear fractional transformation 



CM) 



w*(z) 1 w « (z > = z ~ a ' 



W n [Z) 




- z — a, 

which maps R, U and L onto T, D and E respectively, and has the inverse 

? M-E^M (&«(z)=1-z, 
^a(z) I w* a (z) = a-az. 



Notice that w* a = w a ^, where the ^-involution is now defined by f*(z) = f(~z), 
but nothing similar holds for rh* a . Besides, for the distinguished value a = i, 
C = Ca i s the Cayley transform and ( = ( ao its inverse. 
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Any sequence a = (a n ) n >i in U defines the products (£> n ) n >o as in 
but with the new meaning for ( an . The orthonormalization in L 2 ^ of (-B ra ) n >o 
leads to a sequence (4> n )n>o of ORF with respect to fi with poles in (a n ) n >i, 
which will be called a sequence of ORF on the real line. The study of ORF 
on the real line can be carried out in a completely analogous way to the 
case of the unit circle, so most of the results described for the last ones 
translate directly to the first ones with an obvious change of the meaning in 
the notations. In particular, the sequence (0 n )n>o can be chosen such that 
it satisfies a recurrence like (TTOj) depending on a sequence a = (a n )n>i in D, 
which establishes a surjective application S a : — > D°°, where means now 
the set of probability measures on R. This application is a bijection when 
B n [z) = n^°=i Cn( z ) diverges to zero for z G U, but this is equivalent now to 
Y^=\ I m «n/(1 + \ct n \ 2 ) = oo, which means that the poles can not approach 
too quickly to R. 

Following the same strategy as in the case of the unit circle, we can 
develop a spectral theory for ORF on the real line. The starting point is 
again recurrence ffTUl) written in the form (TT21 . but now the positive factors 
rj a , a G U, are defined by 

(w a (a)\ 1/2 r 

Va=[~^~j =VW 

Both, the expressions for the unit circle and the real line can be combined in 
Va = (zu a (a)/zu ao (a )) 1/2 . 

The form (fl2l) of the recurrence is the key tool to obtain the matrix 
representations with respect to the ORF for the multiplication operator 

/(*) -» 2/(2) 

where \x is the corresponding orthogonality measure on the real line. If 
supp/i is bounded, T M is an everywhere defined self-adjoint operator on L 2 . 
In general, is a densely defined self-adjoint operator on L 2 ^ when the 
function z is finite /i-a.e. (see [2H1 page 259]), that is, when 00 is not a mass 
point of [A. In this case, c p (T M ) = {mass points of /x} and cr(T^) = supp/i 
under the convention that 00 G cr(7]J when T M has an unbounded standard 
spectrum. A way to deal with the case of measures with a mass point at 00 
is to work with the operator multiplication by ( in L 2 , i.e., 

c . 7-2 t2 

/(*)-«*)/(*) 
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This operator is unitary for any measure (« on 1 and verifies the identities 
ffpiS/j) — C( mass points of fi) and = £(supp/i). The matrix represen- 

tations of T M and with respect to the related ORF are related to the 
operator analogs of the new linear fractional transformations ( a . 

To discuss such operator linear fractional transformations it is convenient 
to introduce the notation 

ReT=i(T + T t ), ImT = ^(T-T t ), 

for any densely defined operator T on H. The operator linear fractional 
transformations of interest for ORF on the real line are 

UT) = VA w A (T)- 1 zu* A (T) V 2\ 

( A (T) = r ] - A l w\(T) m A (T)- l VA , 
where 



w A (T) 


= T-A\ 


w* A (T) 


= T-A, 


w A (T) 


= 1-T, 


w\{T) 


= VaAtj^ 1 



VaA^T, 



rj A = Vim A 

and A e B# is such that ImA > e for some positive number e (in short, 
Imi > 0), so that r\ A is bounded with bounded inverse. When A is normal, 
as it is the case related to ORF on the real line, w\ becomes 

zb* A (T) = A - A ] T. 

( A is a bijection of Vh = {T 6 B# : ImT > 0} onto B#, and ( A is its 
inverse. To prove this assertion we start showing that ( A and ( A map U# 
and H>h respectively on B#. The statement for ( A is a consequence of the 
fact that ||T|| < 1 implies w A {T)^ e B#. As for ( A , the result follows from 
the fact that the spectrum of any operator T e B# is included in the closure 
of its numerical range {(x, Tx) : ||x|| = 1}. So, if ImT > 0, then cr(T) C U 
and thus T" 1 e B#. In consequence, tu^T)^ 1 e M H for any T G Uh since 
Im(T- A f ) > ImA > 0. 

On the other hand, using the equality A^rf^A = 2iA+A vf^ A = A f] A 2 A ! , 
we find the identities 

w A (T) r,2\l - (a(T) G(T)t) ifi ViW = 4 ImT, 
^(Tj^^Imair))^ 1 ^^) = 1 -TtT, 
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which prove that (a maps U# on ~B H and (a does the opposite. Moreover, a 
direct calculation shows that, for any T G U# and any S G B#, S = (a{T) 
iff T = (a(S). This completes the proof. 

The above arguments can be easily generalized to see that (a extends to a 
transformation of IJ H = {T G B H : ImT >_0} onto {T G B H : 1 £ a(T)}, ( A 
being its inverse. In consequence, £4 maps Uh-\U# onto {T G T# : 1 ^ 
and does the converse. Furthermore, (T4"6T) also implies that £4 maps the 
set of bounded self-adjoint operators onto the set of unitary operators whose 
spectrum does not contain 1. 

The above properties are verified in particular by the Cayley transform 
for operators, since it is given by ( = Ca=i- Indeed, (a is nothing but the 
composition of the Cayley transform with an operator transformation de- 
pending on A which maps onto theirselves Uh, Vh and the set of self-adjoint 
operators on H. More precisely, taking into account that 

Va^a^Va 1 = Va 1 ( t -Re^VA 1 +h 
r]fw* A {T) Va 1 = Va(T - Re A) rfA 1 - i, 

we obtain 

a(T) = Cte 1 (T-ReA)r ?A 1 ). (47) 

It is known that the Cayley transform extends to a bijection between the 
set of (bounded or unbounded) self-adjoint operators and the set of unitary 
operators whose point spectrum does not contain 1, so the same holds for (a- 
The importance of this property is that it permits us to formulate the spectral 
theory for ORF on the real line, including the case of measures on IR with 
unbounded support since they are associated with unbounded self-adjoint 
multiplication operators. 

Another advantage of relation (f4T|) is that it expresses (a as a product 
of two commutative factors. This provides two equivalent representations of 
(a, namely, 

( A (T) = VA (T - Al)~\T - A) rtf = Va \T - A)(T - A^tja, 

giving rise to two expressions for (a too. From the above result we find that 
Ca^TY = Cai(T^) and (a{TY = ( A i(T^), as in the case of the unit circle. 

Finally, if 3 is a hermitian ideal of B#, similar arguments to those given 
in the proof of Theorem I6. II prove that, for any normal operators A,B G Vh 
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such that AB = BA, the condition A — B G 3 implies the equivalences 

T-Se3 & ( A (T)-( B (S)e3, \/T,SeV H , 

T-Se3 & Ca(T) - Cb(S) G 3, VT, S G Dh, 1 t *(T) U 

The right implication of each case is equivalent to the left implication of the 
other one due to the fact that (a and (a are mutually inverse transformations. 
As we pointed out, when T is unitary the transformation Ca(T) is well defined 
provided that 1 is not an eigenvalue of T. So, the right implication of the 
second equivalence can be formulated in a more general context when the 
operators T, S are unitary. The analogous extension for the right implication 
of the first equivalence, i.e., the case of T, S unbounded self-adjoint operators, 
is not possible because we suppose that 3 is an ideal of B# (as it is the 
case for the classes of perturbations usually considered in operator theory). 
Therefore, we only can assure that 

T-S e3 <= Ca(T) - Cb(S) G 3, VT, S self-adjoint, 

T-S e3 => (a(T) - Cb(S) G 3, VT, S unitary, 1 £ a p (T) U a p (S). 

These results, although weaker than the previous ones, are enough to apply 
perturbative techniques to the spectral theory of ORF on the real line, even 
if the support of the orthogonality measure is unbounded. 

With all these operator tools at hand we can develop the spectral the- 
ory for ORF on the real line following the same steps as in the case of the 
unit circle. In fact, the results for the unit circle are formulated through- 
out the paper in such a way that the translation to the real line is just a 
matter of changing the meaning of the symbols according to the previous 
discussion, together with some other obvious modifications. For instance, 
if fi is a measure on R and (0 n )n>o are the ORF associated with an ar- 
bitrary sequence a. = (a n ) n >i in U, then the sequence (Xn)n>o defined by 
( |32|) . with the new meaning for ( n = ( an , are the ORF associated with 
(a%, &2, «3, «4, • • •)• (Xn)n>o is a basis of L 2 ^ when the odd and even prod- 
ucts B°(z) and B e (z) converge to zero for z G U, but this means now that 

XXa Im a 2 fc-i/(l + \ a 2k-i\ 2 ) = Sfcli Ima 2fc/(l + |«2fc| 2 ) = oo. Also, if 
a = (on)n>i are the parameters of the recurrence for <fi n , the zeros of (f) n 
are the eigenvalues of = C4„(C n ), wnere C — C(°0 and A = A(ot) as 
in the case of the unit circle, but (a.„ is t ne new operator linear fractional 
transformation given in this section. Notice that A n G Uc« and C n G T^n 
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has its eigenvalues in © because they are the zeros of the n-th OP associated 
with the parameters a. Hence, is a well defined matrix of Uc" \ Uc™> 
which agrees with the fact that the zeros of n lie on U. 

Other results for the unit circle can be translated to the real line in a 
similar way, but two of the main results need a special discussion. The first 
one concerns the representation of the self-adjoint multiplication operator 

for a measure fi on K, and the other one is related to the representation 
of the self-adjoint multiplication operator T^v corresponding to the finitely 
supported measure \i v n associated with the PORF Q v n . 

Following the same steps as in Theorem I5.1[ we would find that, if /i 
is a measure on M, for any sequence ex. compactly included in U, the ma- 
trix representation of T M with respect to the ORF (x„)„>o associated with 
(ai, &2, 03, a4, . . .) is U — Ca(C), where A = A(a), C = C(a) and a = S a (fi). 
However, since the matrix C is unitary, we can assure that Ca{C) provides 
a well defined (self-adjoint) operator only when 1 is not an eigenvalue of C. 
That is, in the case of the real line, the matrix representation U = Ca{C) is 
valid provided that 1 <j p {C). To understand the meaning of this condition 
we will relate C to the matrix representation with respect to (Xn)n>o of S^. 
When 1 ^ cr p {C) the matrix of = C(^)t) ^ s C(^0> but, as we will see, an ex- 
pression for the matrix representation of can be obtained for any measure 
fi on the real line, even if it has a mass point at oo. This discussion will lead 
also to a relation between the operator linear fractional transformations in 
the real line and the unit circle. 

Since we are going to consider at the same time the linear fractional 
transformations used on the real line and on the unit circle, in what follows 
we will distinguish between both cases with a superscript K or T respectively. 
Let A G Uh- Due to the properties of the Cayley transform, B = ((A) e D#. 
A direct computation gives 

ImA = (1 - B)~\l - BB^) (1 - Bt)-i. 

Therefore, = \rj^{\ — and, as a consequence of the polar decompo- 

sition, 

r^(l - fi 1 ")" 1 = UrfX, U unitary. (48) 
If we change A by —A* , then B changes to B\ thus, 

nlt(l-B)- x = V7&, V unitary. (49) 
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When A is normal, B is normal too and rf^ = |1 — B\ so U = V' = £b> 
where 

e B = (l-B)|l-B|- 1 . 
In the general case, using (HH1) and (pf9l) . we find that 

Ci(C(T)) = UCl(T)V^, 

llGIlC6 

C(3C0) = C t b(utv1). (50) 

Denoting w = £(z) and 5 = UTV\ a straightforward calculation gives 

- w*{z) Tl = T~~ (<H - **SM 5S) (1 - B)" 1 , (51) 
1 — w V / 

where Tf = (r/^) _1 T r/^ and S|? = (^s) -1 ^ 7 ?!- Since equations f[T5"j) and ([TBI) 
hold for the real line too, the above equality can be written equivalently as 

zwl{Tf) - vjf{Tf) = ^-(wvjUSI) - fitf - B)-\ (52) 

1 — w\ / 

Using (HHD and flU]) we obtain = (1 - B^Tf^l - B). Taking this 
relation into account, a direct computation yields 

1 - Cl(S) = 1 - zul(Sir = {fir^l - T)0€b5(5)- V, 

which implies that 

1 G a(S(5)) ^ 1 G a(T), 1 G 0p(3(5)) & 1 G <x p (T). (53) 

Assume now that a is compactly included in U and \x is a measure on 
R such that 1 ^ cr p (C). From (|50|) we see that the matrix representation 
C(C5(^)) °f ^ can De expressed alternatively as Cb(£sC£b)> with £> = ((A). 
Nevertheless, contrary to C(C5(^))j Cb(^b^^b) is always a well defined (uni- 
tary) matrix, no matter whether 1 is an eigenvalue of C or not, because £g C £g 
is unitary and maps unitary operators into unitary operators. Actually, 
we are going to prove that, if a. compactly included in U, Cb(£bC£b) is the 
matrix representation of with respect to (xn)n>o for any measure \i on R. 
Following similar arguments to those given in the proof of Theorem 15.11 we 
find that, for any measure \i on R, the ORF (xn)n>o satisfy equation ( |35i) too, 
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but substituting C = C\ by C = C*, and zu\, w*J by w\ vj*J- respectively. 
Applying (l5"Tj) and using (ITo]) and (IT5|) we conclude that, for a compactly 
included in U, 

( X o(z) xi(z) •••)(CW-CfeC^)) =0, B = t{A), 

which means that Cb(£bC£b) is the matrix of with respect to (x«)n>o- As 
a consequence of this result and (j551) , we have the equivalences 

1 G a(C) & 1 G ^ 1 G a(^) ^ 1 G C(supp/i), 

1 G crp(C) <^> 1 G 0p(C|(£sC£s)) -v^ 1 G (Tp^) <^> 1 G ((mass points of //). 

Thus, we have reached the following result. 

Theorem 7.1. Let ex. be a sequence compactly included in U, fi a measure 
on R and C = C(a) with a = S a (fi). Then, 

1 G a(C) <=>■ oo G supp/i, 1 G cr p (C) <S> oo is a mass point of ji. 

Therefore, /i is a measure on K iff its related sequence a satisfies 1 ^ cr p (C). 
Thus, U = Ca(C) provides a well defined matrix representation of T M for any 
measure /i on R. Moreover, the measures on R with bounded support are 
characterized by the fact that a is such that 1 ^ &{C). 

In the case of an arbitrary measure \i on R, including the possibility of a 
mass point at oo, we can study the relation /i(a, a) throughout the spectral 
analysis of the matrix representation Cb(£,bC£b) of 5 M or, alternatively, we 
can deal with a pair of operators. More precisely, relation fl52l) implies that 
the spectra of Cb(£bC£b) an d the pair (w^(C),w^(C)) are related by the 
Cayley transform, so 

supp/i = o-(wf(C), w^(C)) = a(C + Ml C\ + A^C Q ). 

Also, the eigenvalues of the pair are the mass points of /x and the eigenvectors 
of the pair with eigenvalue A are spanned by (xo(A) • • • Xn-i(A-) ) VaI/ 2 - 
That is, while the spectral methods that use linear fractional transformations 
of five-diagonal matrices only work for measures on R, their formulation 
in terms of pairs of band matrices are valid for any measure on R. 

Similar results hold too for the finitely supported measures associated 
with the PORF. Given an arbitrary measure /i on R, consider the measure 
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ji v n supported on the zeros of the PORF Q v n = n + t>0*, v G T. As in the case 
of the unit circle, Q v n has n different zeros, but now they lie on R. Besides, 
if u = Can, the matrix representation U^ n ' u ^ = Ctn(^n) °f with respect to 
(Xfc)fc=o * s weu defined whenever 1 ^ cr(C^). Concerning this condition, an 
analogous argument to that of the measure \x proves that 

l e a((%) e> ooe supp/i^, 

i.e., the matrix representation U^ n ' u ^ of T^y is valid for any measure /x^, ex- 
cept for the value v = — 0*(oo)/0 n (oo) which locates a zero of at oo. 
Nevertheless, analogously to the previous discussion, the spectral interpreta- 
tion of the PORF in terms of pairs of band matrices given for the unit circle 
after Theorem 15.91 holds for any PORF on the real line too. 

Concerning the applications of the spectral theory for ORF on the real 
line, from the previous comments we know that, if 3 is an ideal of B^, for 
any sequences ct,f3 compactly included in U and any sequences a, b in D 
such that 1 cr p (C(a)) U a p (C(b)), 

A(a) - A((3), C(a) - C(b) G 3 U(a, a) - U(b, (3) G 3. 

This permits us to extend to ORF on R the applications for ORF on T 
discussed in Section [61 

Equation (1501) provides a connection between the real line and the unit 
circle representations. Let a. = (a n )n>i be a sequence compactly included 
in U, and consider the sequence f3 = (j3 n )n>i in © given by (3 n = ((a n ). 
Following the previous notation we also have a = h so A) = 1- Consider 
two sequences a = (a n ) n >i and b = (6 n ) n >i in © related by 

~2<-2 r"2 t- 1 — fin 



|l-/?n| 



We have the identities C {b) = At£ s C (a) T and C e {b) = T^C e (a) ^gA, where 
B = C(A), A = A(a) and 



™ 1 i ( £062 - ' " 1 °dd n, 

71 1 , 70 = i, 7n = <! ;°; ™ 1 

C1C3 • • • evenn, 

7„_i£„ odd n, 
ln-it; n even n, 



A 

A = I Ax I , A = 1, A n 



(54) 
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Therefore, C(b) = A^ B C(a)^ B A and, thus, equation (loll implies that 

C(W R (a, a)) = A W T (6, (3) A 1 ". (55) 

This relation can be understood taking into account that the ORF on the 
real line and the unit circle are related by the Cayley transform. More pre- 
cisely, 4> n {z) are ORF on the real line iff <f) n (((z)) are ORF on the unit circle. 
If fi is the orthogonality measure on M, the corresponding measure v on T is 
given by v(A) = fi(((A)) for any Borel subset A of T. Also, the parameters 
a n and (3 n associated respectively with the poles of <p n (z) and 4> n {({z)) are 
related by (3 n = ((a n ). Moreover, <p n satisfies the analogue of recurrence ffTUl) 
on the real line with coefficients a n iff (fi n = ■ ■ •^_ 1 ^ n n satisfies such a 
recurrence on the unit circle with coefficients b n = ■ ■ •^_ 1 a n . If Xn and 
Xn are the associated ORF (given by the corresponding version of f !52|) on 
M. and T respectively), then £ n = \ n Xn with \ n as in (|5^|) . Therefore, if a. 
is compactly included in U, the matrix representation U R (a, a.) of with 
respect to (Xn)n>o and the matrix representation U T (b, (3) of T u with respect 
to (Xn)n>o are related by (|53|) . 
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